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ABSTRACT 

We examine the fundamental scaling relations of elliptical galaxies formed through mergers. Using hundreds 
of simulations to judge the impact of progenitor galaxy properties on the properties of merger remnants, we find 
that gas dissipation provides an important contribution to tilt in the Fundamental Plane relation. Dissipationless 
mergers of disks produce remnants that occupy a plane similar to that delineated by the virial relation. As the 
gas content of progenitor disk galaxies is increased, the tilt of the resulting Fundamental Plane relation increases 
and the slope of the R^-Mi, relation steepens. For gas fractions /gas > 30%, the simulated Fundamental 
Plane scalings (R^ oc cr'^^/j"^^) approach those observed in the /T-band (7?e oc cr''^^/J° ™). The dissipationless 
merging of spheroidal galaxies and the re-merging of disk galaxy remnants roughly maintain the tilt of the 
Fundamental Plane occupied by the progenitor ellipticals, approximately independent of the orbital energy 
or angular momentum. Dry merging of spheroidal systems at redshifts z < 1 is then expected to maintain the 
stellar-mass Fundamental Plane relations imprinted by gas-rich merging during the epoch of rapid spheroid and 
supermassive black hole growth at redshifts z « 1 - 3. In our simulations, feedback from supermassive black 
hole growth has only a minor influence on the stellar-mass scaling relations of spheroidal galaxies, but may 
play a role in maintaining the observed Fundamental Plane tilt at optical wavelengths by suppressing residual 
star formation in merger remnants. 

We estimate that « 40- 100% of the Fundamental Plane tilt induced by structural properties, as opposed 
to stellar population effects, owes to trends in the central total-to-stellar mass ratio Mtotai/A^* produced by 
dissipation. Gas cooling allows for an increase in central stellar phase-space density relative to dissipationless 
mergers, thereby decreasing the central Mtotai/A^*- Lower mass systems obtain greater phase-space densities 
than higher mass systems, producing a galaxy mass-dependent central Mtotai/A^* and a corresponding tilt in the 
Fundamental Plane. We account for these trends in the importance of dissipation with galaxy mass in terms of 
the inefficient cooling of coUisionally heated gas in massive halos and dynamically varying gas consumption 
timescales in smaller systems. 

Subject headings: galaxies: formation - galaxies: evolution 



1. INTRODUCTION 

Elliptical galaxies represent a fascinating combination of 
complexity and regularity. A leading theory for the origin 
of e arly-type galaxies is based on mergers of disk galax- 
ies jToomre & T oomr3' 19721 'Toomre"1977|) and likely in- 
volves gas dissipation, star formation, and supermassive black 
hole feedback (Bai-nes 1992; Barnes & Hernauist'199a il996t 
pihos & Hernquist 1994, 1996; Di Matteo et al. 2005) in ad- 
dition to stellar dynamics. Despite their complex origins, 
early-type galaxies obey a regular set of scaling relations that 
connect their photometric and kinematic properties, most no- 
tably the relation between effective radius R^, central stellar 
velocity dispersion cr, and average central surface brightness 
le known as the Fundamental Plane (FP: lDressler et al]l 19871: 
iDiorgovski & Davisbl987.) 



iFaber et al."1987. hereafter the "virial scaling"). 

While the observational determination of the FP was orig- 
inally motivated as a precise distance indicator to improve 
upon the previously known luminosity (L) - velocity dis- 
persion (cr) relation (Faber & Jackson 1976), the importance 
of the FP scalings and its correspondingly small scatter for 
theories of elliptical galaxy formation was also realized. 
The first observationally determined FP scaUngs (a ^ ^ 
1.4,/? - 0.8-0.9, at optical wavelengths: lDressler et al.ll987t 
Diorgovski & Davis 1987) differed from the virial scaling, in- 
dicating a "tilt" relative to the expectation for homologous 
systems. The FP tilt implied that the mass-to-light ratio M/L 
likely varies as a function of galaxy mass or luminosity as 



(2) 



R^ CX cr"/: 



(1) 



The virial theorem can be used to calculate this relation 
for homologous systems, which gives a = 2, /5 = 1 (e.g., 
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within the range 7 « 1 /5 - 1 /4. iFaber et all Jl987ft noted that 
deviations from the FP can be induced hy M/L variance ow- 
ing to e.g. metallicity or age trends in stellar populations, 
dynamical or structural properties, and the relative distribu- 
tion of dark and baryonic matter In principle, each of these 
effects may also introduce a systematic tilt into the observed 
FP if they vary as a function of elliptical galaxy mass. 

The purpose of the current paper is to gauge the importance 
of various contributions to the tilt in the observed FP in the 
context of the scenario where elliptical galaxies form from 
mergers. Using hundreds of simulated galaxy mergers that 
include the physics of gas cooling, star formation, supernova 
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feedback, and black hole accretion and feedback, we deter- 
mine that gas dissipation may significantly contribute to the 
tilt of the observed FP, in addition to tilt induced by M/L 
trends from stellar populations. We propose that elliptical 
galaxies initially form in gas-rich mergers from disk galaxy 
progenitors whose gas fractions exceed /gas ~ 30%, and show 
that these remnants display substantial FP tilt. 

We connect the origin of this tilt to the central stellar phase- 
space density of the remnants. In small mass systems where 
dissipation is most important, the stellar phase-space density 
of remnants increases substantially in high-gas fraction merg- 
ers. The central stellar phase-space density in massive ellip- 
ticals remains similar in mergers with varying gas fractions, 
with their stars on average obtaining lower phase-space den- 
sities than reached in smaller systems. This mass-dependent 
phase-space density trend translates into a mass-dependent 
trend in the ratio of total mass to stellar mass Mtotai/-^* in 
the central regions of ellipticals and a corresponding tilt in 
the FP. We then explain this mass-dependent importance of 
dissipation in terms of the inefficient cooling of collisionally 
heated gas in massive ellipticals and the dynamically vary- 
ing gas consumption timescale of smaller systems. It is inter- 
esting, and possibly significant, that the gas fraction required 
to reproduce the observed tilt in the FP is, as discussed by 
iHernquis? ("1993) similar to that needed for mergers of disks 
to yield the high central phase space densities of ellipticals. 

In accord with other work (e.g., Cacelato et al. 1993 
i Dantas et alJl2003t iNipoti et al J 12003: .Bovlan-Kolchin et all 
l2005h . we find that subsequent dissipationless merging be- 
tween spheroidal galaxies roughly maintains the FP tilt. 
Moreover, we find that a single generation of dissipationless 
re-merging of remnants will induce scatter in the Mbh-c rela- 
tion but will not destroy the correlation, as spheroidal galaxy 
mergers also do not dramatically alter the R^-M^, relation 
when viewed as a mass-sequence. Possible dry merging be- 
tween spheroidal galaxies at redshifts z < 1 as indicated by re- 
cent observations (Bell et al. 2005; van Dokkum 2005) is thus 
not expected to destroy tight FP orM^n-a relations generated 
during spheroid formation through gas-rich mergers at higher 
redshifts. 

This paper is organized as follows. We review the observa- 
tional and theoretical work on elliptical galaxy formation and 
scaling relations in ^ We present our methodology in |3l™d 
our results in |3 We discuss the implications of our work in 
Inland summarize and conclude in ^ Throughout, we adopt 
a flat ACDM universe with = 0.3, = 0.7, Ob = 0.04, and 
a Hubble parameter Hq = lOO/i km s"' Mpc"' with h = 0.7. 

2. REVIEW OF ELLIPTICAL GALAXY SCALING LAWS 

Scaling laws describing the regularity of the properties o f 
elliptical galaxies have been known since at least lFish] J 19641) . 
who reported a relation between their potential energy W and 

mass M^, as W oc M^^. Even earlier, Ide VaucouleursI J 19481 
hereafter dV) had demonstrated that ellipticals generally have 
a surface brightness profile log/(r) oc r'/'*, at least over some 
range in radius r, and implications of the dV profile for the 
mass-dependent properties of ellipticals were recognized be- 
fore the Fish ( 1964) paper (e.g., Poveda,1958J . 

Sandage ( 1972) found a color-magnit ude relation for Virgo 
and Coma cluster ellipticals (see also 'Stebbins & WhitfordI 
[^52; de Vaucouleurs 1961). Faber (1973) discovered a sim- 
ilar color-magnitude trend in Local Group and cluster ellip- 
ticals, as well as an absorption-line strength-magnitude rela- 
tion, and suggested that these properties of elliptical galax- 



ies depend primarily on galaxy luminosity. iFaber & JacksonI 

( 1976) established the relation between luminosity L and ve- 
locity dispersion cr, providing further evidence that ellipticals 
follow a regular sequence as a function of mass. Kormend^ 

( 1977) showed that the surface brightnesses and effective radii 
of ellipticals correlate with gal axy lu minosity and with one 
another (see also Binggeli et al .119841) . 

Important early indications that a second parameter (in ad- 
dition to mass) governs the pr opert ies of elliptical s came 
with the Terlev ich et al] ('1981) and iTonrv & DavisI (119811) 
work that implied a correlation of L - tr and absorption-line 
strength- luminosity relation (MgT- L) residuals. Wh ile these 
findings were later contradicted (e.g. lDressledll984 l). that el- 
liptical galaxies were not a one-parameter family remained an 
important possibility. 

The discovery of the FP (iDressler et alJ Il987t 
Diorgovski & Davis 1987) definitively revealed that el- 
liptical galaxy properties are set by at least two parameters. 
Specifically, ellipticals were found to obey a relation between 
^e, cr, and 4 as given by Equation dD, with less than half 
the scatter of the Faber & Jackson ( 1976) L-a relation. The 
small scatter of the FP was immediately noticed, implying 
that the process of elliptical galaxy formation must result in 
a very regular mass-sequence. Faber et al. ( 1987) noted that 
mass-to-light ratio {M/L) variations can influence observa- 
tions of the FP by inducing tilt relative to the plane defined 
by the virial relation. Much of the subsequent work on the 
FP has centered around possible causes of M/L variation or 
other origins for tilting the FP relative to the virial plane. 

Numerous subsequent observations verified and im- 
proved the FP relation for ellipticals ( Lucev et afl 
1991a b; de Carvalho & Diorgovski 1992; Bender etaD 

1992; J0rgensen et al. 1992; Guzman et al^ ,d^2S 

J0rgensen et al. 1993; SagUa et al. 1993; Bender etai.lll994t 
Prugniel & Simien 1994; Pahre et al. 1995; J0rgensen et alJ 
1996; Prugniel & Simien 1996; Busarello et al. 1991 
Graham & Colless 1997; Bender etal. 1998; Pahre et an 
|1998b.a£ Mobasher et aL .1999; Krona witter et all 
|2000t iGerhard et alJ 12001 iBern ardi et a l.1 I2003at 
'Padmanabhan et al."2004'; 'Woo et al."2004; Cappellari et aiJ 
2005 ) and extended the FP determinations to higher redshifts 
(Franx 1993; van Dokkum & Franx 1996; Kelson etaj 
1997; Schade etal. .1997: van Dokkum et al.. .199^ 
.TOrgensenetalJ 1l999t iTreuetalJ l1999t iKelsolT^ 
2000a b"ct 'Kochanek et al.' '200ff, 'Tr eu et al.1 [2001; 
van Dokkum & Franx 2001 ; van Dokku m et alJ |2001 ; 
llreuetal. 2002; Gebhardt et al. 2003; Rusin et al. 200:^; 
Ivan Dokkum & Stanford. ,2003: .van Dokkum & Ellis 2003 i 
van der Wei et al. 200^7 While the specific details may vary, 
these works generally conclude that 

• A tight fundamental plane between the elliptical galaxy 
properties 7?e, cr, 4 exists, and extends in some form to 
at least redshift z ~ 1 . 

• Ellipticals are old, with formation redshifts z > 1, and 
their color evolution, which controls the FP normal- 
ization, is roughly consistent with passive evolution of 
their stellar populations. 

• Some of the FP tilt must originate from the change in 
stellar population M/L with galaxy mass, but the extent 
to which stellar populations contribute to M/L-induced 
FP tilt is debated. 



Gas Dissipation & the FP 



3 



Of specific interest to modelers are the observations that made 
definite statements about the influence of structural or kine- 
matic nonhomology on M/L and the FP tilt, especially those 
that conclude directly that these nonho mologies are either 
unimportant or of minor significance (Gerhard et al. 2001; 
l^^ppellari et al. 2005), significant (e.g. Padmana bhan et al„ 
120041) . or of possibl e but as y et not fully determined impor- 
tance (e.g.l Pahre et al.ll998bh . 

Moreover, elliptical galaxy properties related to the 
FP or its projections have been extensively observed. 
These include studi es of the photom etric profiles (e.g. 

Vaucouleurs 1948| [Ser^i(|:||ip^ ^ iKorm endv & Illingworth 
[1982^ Burkert 1993t ICaon et alJ Il993i). internal" kine- 
matic structure (e.^ . iBhrnevT 119781 iDav iesetal. 1983 
iDavies & Birkinshaw 1988; Bender etal. 1994), metallic- 
ity (deCarvalho&Diorgovski 1992; Bender_et al. 1991 
[1996; Bernai-di et al. l f]"998l l2603bl). and the Mrh-ct rela- 
tion (e.g. Gebhardt et all 1200 0^. 'Ferrarese & Merritll 120001 
[Tremaine et al. 2002|1 An important, related prop erty is 
the size-stellar mass relation (Shen et al. 2003; Truiill oet al.l 
12004; Truiillo & Aguerri 2004; Mcintosh et aL^QQSj, that in- 
dicates a power-law correlation between some characteristic 
galaxy size and the stellar mass or luminosity. 

Interpreting these observations has been the focus of 
various theoretical efforts (see, e.g. iBarnes & HernauistI 
[1992, for a description of early results). Notably, sim- 
ulations of the formation of ellipticals and their prop- 

the merger hypothesis have 
variety of p r ogenitor mod- 
(White 197?; Capelato et all 



erties in the context of 
been performed using a 
els including spheroidal 
1.995; Dantasetal. 2003; 
2003; NiDotietal. 2003), 



Gonzalez-Garcia & van Albadi 
disk ( Toomre & Toomre 1972; 
Gerhard.l98L Farouki & ShapirQ.1982; .Faro uki et al. 1983 
Bames&IIernc|ui p I199U iBarnesI 119921 iHernauisti 1199 



]^93VHernauist et al."1993';'Mihos & Hernauist 1994, 1996 ; 
Hibbard & van Gorkom 1996; Bekki 1998; Dubinski 199| 
Naab et al. 1999; Naab & Burkert 2001; Aceves & Velazaued 
2005.). and cosmological systems (e.g. Aarseth & Fall 1983 
Saiz et all [2 0041 ). Analytical models of elliptical galaxies 
have also been formulated (Hernauist 1990; Ciotti .19911 



[Ciotti et al. 1996; Ciotti 1996; Ciotti & Bertin 1999), aiding 
the interpretation of both the observational and simulation 
results. In what follows, we combine features of many of 
these previous theoretical endeavors by simulating mergers 
between dissipational and dissipationless disk galaxies, 
spheroidal systems, and merger remnants to determine their 
fundamental scaUng relations. 

3. METHODOLOGY 

Galaxy merging combines much complex physics, includ- 
ing e.g. the collisionless dynamics of dark matter and stars, 
gas dissipation, star formation, and feedback from supernovae 
and black hole growth. While each of these processes may in 
principle be essential for determining merger remnant proper- 
ties, their relative importance has not been fully established. 
By comparing the structure of merger remnants in simulations 
which systematically include or exclude various processes, 
we attempt to both test the merger hypothesis and identify 
the most important physical mechanisms. 

To this end, we perform a set of hundreds of simulations of 
galaxy mergers. Our suite consists of three categories of disk- 
disk mergers and one class of spheroid-spheroid mergers. The 
first category of simulations, referred to as "dissipationless," 
includes only simulations of disks consisting of collisionless 



stars and dark matter The second category, termed "dissipa- 
tional" mergers, consists of simulations that account for gas 
cooling, star formation, the physics of the interstellar medium, 
and supernova feedback. We refer to the third category of cal- 
culations as "full-model" simulations which, in addition to the 
processes accounted for in the "dissipational" mergers, also 
include supermassive black hole growth and feedback. All 
three categories of disk-disk mergers are derived from a stan- 
dard set of galaxy models to enable a direct comparison of 
physical processes. 

To perform our numerical simulations, we utilize the 
GADGET2 code (Springel 2005). GADGET2 use s flie 
smoothed particle hydrodynamics (SPH) formalism JLucvl 
[l977; Gingold & Monaghan 1977) in its entropy-conserving 
formulation ( Springel & Hernauist 2002) to calculate the dy- 
namical evolution of the gas and a tree-bas ed method to 
comp ute gravitational forces between particles JBarnes &Hutl 
fT986l) . 

All pro genitor galaxies were created using the methods de- 
scribed in Springel et al. (2005), which allow for the genera- 
tion of stable equilibrium galaxy models. Each galaxy con- 
tains an extended dark matter halo, and may also consist of 
a stellar disk, gaseous disk, stellar bulge, and a supermas- 
sive black hole particle. The collisionless components of 
the galaxy models are required to satisfy the Jeans equations 
while the structure of the gas component is determined by the 
equation of hydrostatic equilibrium and an integral constraint 
on the surface mass density. 

3.1. Dissipationless Disk Simulations 

The dissipationless disk progenitors consist of an exponen- 
tial stellar disk embedded in a dark matter halo with virial 
velocities in the range Vvir = 80 - 500 km s~'. We initial- 
ize the size of the disk according to the iMo etal J figgl 
formalism for di ssipation al disk galaxy forma tion (see also 
iFaU & Efstathioutl980;.Blumenthal et alJl986h assuming the 
disk contains a fraction of the total galaxy angular mo- 
mentum equal to its mass fraction, which we set to = 
Mdisk/Mvir = 0.041 to match the Milky Way-like model used 
by Springel et al. (2005). The disk scalelength is then de- 
termined by the galaxy spin A and the dark matter halo con- 
centration Cvir. We adopt A = 0.033, which is near the mode 
of the redshift- and mass- independent distribution of dark 
matter halo spins measured in cosmological N-body simu- 
lations (Vitvitskaet al. 2002). For the Navarro-Frenk- White 
halo concentration Cvh (Navarro et al. 1997, NEW), we adopt 
the mass- and redshift-dependent dark matter halo concentra- 
tions measured bv lBuUock et al.l (1200 10 



Cvir(Mvi,-,z)«9 



\A^coll,0 



-0.13 



(3) 



where Mcou.o '~ 8 x lO'^/z 'M© is the linear collapse mass at 
redshift z = 0. In all cases, the stellar disk scaleheight hi = 
0.2rd, similar- to the Milky Way (c.f. Siegel et al. 2002). We 
model the dark matter halo with a .Hernauist. 11990.) density 
profile of the form 



Ph(r) = ■ 



(4) 



27r r{r+af' 

where the scalelength fl(Cvh ) maps the 'Hernauist' ('1990") pro- 
file param eters to the appropria te NEW halo parameters (for 
details, see lSpringel et al.l2005ft . 
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Table 1. Galaxy Mergers 



Model Progenitor Redshift Gas Fraction ISM Pressurization Pericentric Separation # of Simulations 

Vvir[kms"'] Z /gas C/EOS '-peri 



"Full-Model" Simulations with Black Holes 



Local 


80, 115, 160, 226, 320, 500 





0.4, 0.8 


0.25, 1.0 


2rd 


24 


Intermediate— z 


80, 115, 160, 226, 320, 500 


2,3 


0.4, 0.8 


0.25, 1.0 


2rd 


48 


High-z 


115, 160, 226, 320, 500 


6 


0.4, 0.8 


0.25, 1.0 


2rd 


20 


Halo Concentrations 


160 





0.4 


1.0 




5 


Disk Orientation 


160 





0.4 


1.0 


Tabled 


14 


Orbital Configuration 


160 





0.4 


1.0 


TablelH 


18 


Dissipational Simulations 


Local 


80, 115, 160, 226, 320, 500 





0.4, 0.8 


0.25, 1.0 


2rd 


24 


Intermediate— z 


80, 115, 160, 226, 320, 500 


2,3 


0.4, 0.8 


0.25, 1.0 


2'-d 


48 


High-z 


115, 160, 226, 320, 500 


6 


0.4, 0.8 


0.25, 1.0 


2rd 


20 


/gas Runs 


80, 115, 160, 226, 320, 500 





0.01,0.025,0.05 


0.25, 1.0 


2rd 


72 








0.1,0.2, 0.4 








Dissipationless Simulations 


Local 


80, 115, 160, 226, 320, 500 





0.0 




2ri 


6 


Intermediate— z 


80, 115, 160, 226, 320, 500 


2,3 


0.0 




2ri 


12 


High-z 


80, 115, 160, 226, 320, 500 


6 


0.0 




2rd 


6 


Wide Orbit 


80, 115, 160, 226, 320, 500 





0.0 




OAR,,, 


6 


Wide Orbit, Int.-z 


80, 115, 160, 226, 320, 500 


2,3 


0.0 




OAR,,, 


12 


Wide Orbit, High-z 


80, 115, 160, 226, 320, 500 


6 


0.0 




OAR,,, 


6 


Bulge, Local 


80, 115, 160, 226, 320, 500 





0.0 




2rd 


6 


Bulge, Intermediate— z 


80, 115, 160, 226, 320, 500 


2,3 


0.0 




2rd 


12 


Bulge, High-z 


80, 115, 160, 226, 320, 500 


6 


0.0 




2rd 


6 


High-Res 


80, 1 15, 160, 226, 320, 500 





0.0 




2rd 


6 


Spheroid Simulations 


Local 


80, 115. 160, 226, 320, 500 





0.0 




0.025Rvir 


6 


Wide/Elliptical Orbit 


80, 115, 160, 226, 320, 500 





0.0 




0.4Rvi, 


6 


Full-Model Re-mergers 


80, 115, 160, 226 





0.4 


0.25 


0.05Rv,r 
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Following'R obertson et al.N2005lJ) . we scale the progenitor 
galaxy properties to approximate the structure of disk galax- 
ies appropriate for redshifts z = 0, 2, 3, and 6. Varying the 
progenitor galaxies in this manner enables us to determine the 
impact of redshift-dependent galaxy properties on the scaling 
laws of ellipticals. Keeping the virial velocity Vvir fixed with 
redshift, we scale the progenitor virial mass and virial radius 
using the relations 



lOGH(z) 



Vv, 



lOHizy 



(5) 



(6) 



where H{z) is the Hubble parameter To suitably resolve the 
forces between particles in models of higher-redshift systems, 
we reduce the gravitational smoothing by (1 -fz)"'. The halo 
concentrations also vary with redshift and mass according to 
Equation (|3} and the disk scalelengths decrease with redshift 
through their dependence on 7?vir and Qir- We note that the 
redshift-dependence of disk scalelengths agrees well with the 
distribution of disk scalelengths seen out to redshift z « 1 (e.g. 

iRavindranath et al..2004: Harden et al 2005) . 

At each redshift, we consider three separate types of dis- 
sipationless disk mergers. First, we examine equal mass 
mergers of pure disk galaxies on prograde-prograde copla- 
nar parabolic orbits with the pericentric passage distance set 
to Tperi = 2rd. The galaxies each contain 60,000 dark mat- 
ter and 80,000 stellar disk particles. Second, we repeat each 



merger with a wider parabolic orbit increased to rpeii = 0.47? vir 
to judge the effect of increased orbital angular momentum on 
the dissipationless merger remnants. In addition, we repeat 
each nearly-radial merger with bulge components included in 
the galaxies with mass fraction m\, = Mbuige/A^vir = .1367 
to mat ch the Milky Way-like mod el used in Soringe l et alJ 
J2005I) . We model the bulges with a iHernauisll (k 1 99 ) density 
profile form (see Equation]^ where we set the bulge scale- 
length b = 0.2rd. Each bulge contains 20,000 particles, with 
the number of disk particles reduced to 60,000 to maintain 
the same mass resolution. Finally, we re-run all the pure disk 
simulations at redshift z = with higher resolution dissipation- 
less models with 180,000 dark matter particles and 120,000 
disk particles to examine issues related to numerical resolu- 
tion. While we discuss these tests in more detail in ^ we 
note here that the large set of dissipationless simulations pro- 
duces results very consistent with the restricted set of higher 
resolution runs. In all we perform 78 dissipationless simula- 
tions, and we provide a complete listing in TableQ] 

3.2. Dissipational Disk Simulations 

To gauge the impact of dissipational gas physics on the 
properties of merger remnants, we perform a suite of disk 
galaxy mergers that include gas cooling, star formation, and 
supernova feedback. The dissipational disk progenitors con- 
tain e xponential gaseous and stellar disks and Hernaui^ 
( 1990) dark matter halos, with their disk sizes determined by 
the Mo et al. (1998) formalism as described in ^3.11 The ver- 
tical structure of the gaseous disks are determined by an inte- 
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Table 2. Orbital Variations 



Models 


di 




02 




' peri 




[deg] 


[deg] 


[deg] 


[deg] 


[/r' kpc] 




Disk Orientations 


b 


180 











5.0 


c 


180 





180 





5.0 


d 


90 











5.0 


e 


30 


60 


-30 


45 


5.0 


f 


60 


60 


150 





5.0 


g 


150 





-30 


45 


5.0 


h 














5.0 


i 








71 


30 


5.0 


j 


-109 


90 


71 


90 


5.0 


k 


-109 


-30 


71 


-30 


5.0 


I 


-109 


30 


180 





5.0 


m 








71 


90 


5.0 


n 


-109 


-30 


71 


30 


5.0 


o 


-109 


30 


71 


-30 


5.0 


P 


-109 


90 


180 





5.0 


Orbital Configurations 


el 


30 


60 


-30 


45 


2.5 


el 


30 


60 


-30 


45 


10.0 


el, 


30 


60 


-30 


45 


15.0 


e4 


30 


60 


-30 


45 


20.0 


e5 


30 


60 


-30 


45 


40.0 


e6 


30 


60 


-30 


45 


30.0 


hi 














2.5 


hi 














10.0 


h3 














15.0 


h4 














20.0 


h5 














40.0 


h6 














30.0 


kl 


-109 


-30 


71 


-30 


2.5 


kl 


-109 


-30 


71 


-30 


10.0 


k3 


-109 


-30 


71 


-30 


15.0 


k4 


-109 
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gral constraint from the surface mass density and the require- 
ment of hydrostatic equilibrium within the galaxy potential. 

The thermal properties of the gas are determined us- 
ing the mu ltiphase interstellar medium (ISM) model of 
ISpringel'&^Hern auis? (120031). Star formation is prescribed 
in the manner of Springel & Hernquist (2003), constrained 
to approximate the Schmidt ( 1959) law for disk galax- 
ies as measured by Kennicutt (1998), including a density 
threshold. Below this threshold the gas is modeled as a 
single-phase medium which is not star-forming. Dense gas 
above the threshold is modeled as a hybrid of cold, dense 
clouds embedded in a hot, diffuse medium as envisioned by 
|McKee & Ostriker ( 1977). The temperature of the hot phase 
is set by supernova feedback and the efficiency of cloud evap- 
oration, and has an energy per unit mass that far exceeds the 
cold phase. Even though most of the gas by mass is cold, the 
high temperature of the hot phase more than compensates for 
its small mass fraction, acting to pressurize the star-forming 
gas, and leading to an effective equation of state Ppffip) that is 
stiffer than isothermal (for a numerical fit, see Robertson et alj 
l2004h . The multiphase model of Scrinael & Hernauist (2003) 
has been generalized by Spri ngel et al.. (.2005 ) to allow for 
an effective equation of state parameter qEos that linearly 
interpolates between an isothermal gas (^eos = 0) and the 
fully-pressurized multiphase ISM model (q'eos = !)■ In- 
creasing ^Eos improves the dynamical stability of the gas 



and can prevent 'Toomrd (119641) instability even in gas-rich 
system s ( Springel & Hernauist 2003; Robertson et al J 120041: 
ISpringel & Hernauist.2005i: .Robertson et alj2005aiK 

For our dissipational models, we re-run theVvir = 80-500 
km s"' pure disk merger simulations from ^3. II with two gas 
fractions of /gas = 0.4,0.8, each with two equation of state pa- 
rameters ^Eos = 0.25, 1.0, at redshifts z = 0, 2, 3, and 6. Each 
progenitor galaxy has 60,000 dark matter particles, 40,000 
stellar disk particles, and 40,000 gas particles. The systems 
are merged on prograde-prograde parabolic coplanar orbits 
with Tperi = 2rd. For an equation of state parameter qEos = 0.25 
we also systematically vary the gas fraction of z = progen- 
itors using /gas = 0.01, 0.025, 0.05, 0.1, 0.2 and 0.4. Our 
dissipational simulation category has a total of 164 runs, with 
the complete list of simulations provided in Tabled 

3.3. Full-Model Disk Simulations 

Our full-model category simulati ons include the complete 
physical model presented in Sprin gel et alJ ( l2005h . account- 
ing for gas cooling, star formation, supernova feedback, 
t he S pringel & Hernauist (2003) ISM model described in 
^3.21 and a prescription for supermassive black hole growth 
and feedback. The supermassive black holes are included 
as "sink" particles, with seed masses of 10^ HT^Mq. The 
black holes are allowed to grow according to spherical 
Bondi-Hoyle-Lyttleton accretion ( Hovle & Lvttleton 19391 
fBondi & Ho vle 1944; Bondi 1952). The mass accretion rate 
M is determined from the density and sound speed of the 
gas near the black hole. A fraction ef = 0. 1 of the accretion 
rate is radiatively released, of which a fraction 77theim = 0.05 
is coupled as thermal feedback into gas within an SPH ker- 
nel smoothing length of the black hole. The strength of the 
thermal coupling is comparable to the thermal feedback cou- 
pling of supernova en ergy used in cosmological simulations 
(e.g. Abadi et al. "2003), and reproduces the Mbr-ct relation 
observed locally (Di Matteo et al. 2005). Using the same full- 
model simulations presented here, Ro bertson et a l. (2005lJ 
have also demonstrated that this model for black hole accre- 
tion and feedback should preserve the power-law scaling of 
the Mbh-o" relation between redshifts z = 0-6. 

The full-model simulations augment the pure disk dissipa- 
tional merger simulations from ^3. 21 with supermassive black 
hole growth as described above. Each progenitor galaxy has 
40,000 stellar disk particles and 40,000 gas particles, and are 
merged on prograde-prograde parabolic coplanar orbits with 
'"peri = 2rc|. The models are calculated for Vyir = 80-500 km 
s"' galaxies at z = 0, 2, and 3 and Vvir = 115-500 km s"^ 
galaxies at z = 6. Simulations are performed with gas frac- 
tions of /gas = 0.4,0.8 each with equation of state parameters 
of <?EOS = 0.25, 1.0. Furthermore, we run 14 variations of the 
Vvir =160 km s"' progenitor merger where we change the disk 
orientation according to the method of Barnes ( 1992) to char- 
acterize the effects of disk alignment (see Table|2li, and run an 
additional 18 simulations where for 3 different orientations 
we vary the pericentric passage distance. We also run a set of 
5 additional simulations of the Vvh = 160 km s"' halo where 
we simulate dark matter halo concentrations of Cvir = 5, 7, 9, 
11, and 13. In all, we perform a total of 129 full-model sim- 
ulations and a complete list of these runs is provided in Table 

3.4. Spheroid Simulations 

To explore the properties of remnants formed by the merg- 
ing of spheroidal systems, we also perform a suite of equal 
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mass spheroid-spheroid mergers. The spheroid progenitors 
consist of st ellar spheroids em bedded in dark matter halos, 
both having iHernauistI ( IT990l) profiles. We assume a stellar 
mass fraction of = 0.05 and the concentrations of the halos 
are adjusted to account for the mass-dependence measured 
in cosmological simulations (see Section ISjTl. The s i zes of 
the stellar spheroids are set to follow the IShen et alJ J2003h 
Rg-Mi, relation for massive galaxies as 



R^ = A.l6 



lOiiM, 







kpc, 



(7) 



(see also iBovlan-Kolchin et aljf2005h . We vary the circular 
velocity of the halos between Vvii = 80 - 500 km s"' . The 
spheroid galaxies merge on either nearly-radial parabolic or- 
bits with Tperi = 0.025^vir or wide elliptical orbits with ellip- 
ticity e = 0.5 and rpeii = 0.47?vir- Each spheroid progenitor has 
1,200,000 dark matter and 80,000 stellar particles. A com- 
plete list of the spheroid-spheroid mergers is provided in Ta- 
ble[T] 

To determine the impact subsequent re-merging between el- 
lipticals formed from disk galaxy mergers might have on the 
structural properties of the remnants, we also re-merge disk 
galaxy remnants from the z = full-model simulations. These 
remnants are merged on parabolic orbits with rperi = 0.05/?vir, 
where 7?vii is the virial radius of the original disk progenitor. 
We re-merge remnants with progenitor galaxy circular veloc- 
ities in the range Vyir = 80-226 km s~' . A complete list of the 
re-merger simulations is provided in Tabled 

3.5. Analysis 

Each simulation is evolved until the merger is complete and 
the remnants are fully relaxed, requiring integrations of typi- 
cally 2-4 Gyr The remnants are then kinematically analyzed 
by measuring the half-mass stellar effective radii R^, the av- 
erage one- dimensional velocity dispersion a within a circu- 
lar aperture of radius R^, and the average stellar surface mass 
density 4 measured within R^ as /g = M*(r < Rf.)/TTRl. The 
quantities /?e, cr, and 4 are averaged over 100 random sight 
lines to the remnant. 

Once the FP parameters 7?e, cr, a nd 4 are determined, we 
employ the direct fitting method of Bernardi et al. ( 2003a) to 
determine the best-fit FP scalings. The direct fitting method 
seeks to minimize 



A = log/?e-alogcr-/31og/e-5. 



(8) 



where a and f3 are the FP scaling indices defined by Equation 
{I). The minimization of A requires 



(c^iVlv- 



■c^irO-iv) 



vv 



crLcr?v) 



-^vv 



5=(l0g/?e)-«(l0gfT)-/3(l0g/e), 



(9) 



(10) 



(11) 

where the average (logX) over the data samples considered 
is defined as 



{\ogX)=Y,^QgXi/N 



(12) 



and the co-variant dispersion is defined as 

(logX,- (logX)) (logi^.- (logy)) /N. (13) 



Tables. Best-Fit Scalings 



Models 




Fundamental Plane 








a 




A 






Dissipationless 


2.00 


1.01 


l.OOitO.Ol 


0.018 


0.45 ±0.03 


Dissipational 


1.58 


0.80 


0.80 ±0.01 


0.065 


0.57 ±0.02 


Full-Model 


1.55 


0.82 


0.79 ±0.01 


0.062 


0.57 ±0.02 


Gas Fraction /gas Runs 


Aas = 0.01 


1.81 


0.75 


0.97 ±0.01 


0.009 


0.41 ±0.01 


/gas = 0.025 


2.11 


0.74 


0.96 ±0.01 


0.011 


0.41 ±0.01 


/gas = 0.05 


2.07 


0.64 


0.95 ±0.01 


0.011 


0.41 ±0.01 


/gas = 0.1 


2.01 


0.61 


0.92 ±0.01 


0.014 


0.42 ±0.01 


/gas = 0.2 


1.89 


1.20 


0.89 ±0.02 


0.024 


0.44 ±0.01 


/gas = 0.4 


1.64 


1.07 


0.83 ±0.02 


0.033 


0.51 ±0.03 



The mean-squared scatter about the direct best-fit plane can 
then be characterized by the quantity 



2 2 2 
''ll'^RR^VV" 



(14) 



4 ,^2 2 2 



X [<7„a- 



II" vv" 



^iv; 



When appropriate we compare the quantity (A^) with the 
scatter determined from observational samples. 

The remnant properties are also compared to the R^ -M* re- 
lation, which has been measured observationally in the SDSS 
JShenetalJl200l and may be represented by the power-law 
form 

/JeOcMi". (15) 

The 7?e -Mi, relation allows for a useful comparison of rem- 
nant sizes for simulations with differing angular momenta, 
progenitor redshifts, ISM physics, or gas dissipation. In addi- 
tion, we find that the relative location of remnants in the FP 
can often be related to the impact of different physical pro- 
cesses on the effective radii. The /?e-M* relation also serves 
as a useful calibration for our method to measu re the FP prop- 
erties of the remnants. As mentioned in ^3.41 a subset of our 
simulation suite involves the merging of equilibrium models 
of spheroids initialized to satisfy the Shenet al. (2003) rela- 
tion. As discussed in further detail in ^4.41 the analysis tech- 
nique used to measure the FP properties of remnants accu- 
rately recovers the Shen et al. (2003) relation when applied to 
the spheroid model progenitors and affirms our ability to de- 
termine simulated remnant properties with reasonable fidelity. 

Typically the root-mean-squared variation in the velocity 
dispersion and effective radius along different lines-of-sight to 
a given remnant are Act k, 20 km s"' and A7?e ~ 0.5/i~' kpc, 
respectively, for remnants formed through mergers of Vvir = 
160 km s"' progenitors. The magnitude of the line-of-sight 
variation scales with the velocity dispersion or effective ra- 
dius, but the fractional line-of-sight variation remains similar 
with Act « 0.1 (cr) and AT^e « 0.1 (7?e) roughly independent 
of progenitor mass. The combined quantity cr^/j' has a typi- 
cal fractional line-of-sight variation of Acr^/j' w 0.3 (cr^/j') 
roughly independent of progenitor mass and under this varia- 
tion galaxies tend to move along relatively short trajectories in 
the /^ec/e-space that do not deviate far from the mean FP. Cer- 
tain orbits, such as the head-on orbit considered for spheroidal 
mergers in ^4.41 can induce larger line-of-sight variations but 
the average variations reported above are much more frequent 
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in our remnants. The impact of these relatively small errors on 
estimates of the power-law correlations of the FP and R^-Mi, 
relations are minor and comparable to deviations induced by 
either orbital variations or the structural properties of progen- 
itors, and therefore additional scatter induced by the line-of- 
sight variations of individual systems is not considered by our 
subsequent analysis. However, when plotting either the FP or 
the R^-Mi, relations, we provide sample error bars to show 
the rough deviations any single galaxy could make from their 
mean positions due to line-of-sight variations in cr or /?e. 

3.6. Comparison with Observations 

We adopt the approach of plotting the FP properties of sim- 
ulated remnants in the R^-a^I'^ virial plane coordinate sys- 
tem. The tilt of the FP relative to the virial plane can be quan- 
tified through the power-law relation 

R,^{a^i;')\ (16) 

where A = 1 indicates an alignment of the FP with the virial 
plane. When appropriate, we measure this estimate of the FP 
tilt relative to the virial plane by linear least-squares fitting. 
Our choice of coordinate systems is not unique, and alterna- 
tive rep resentations of th e FP include the K-space coordinate 
system ( Bender et al. 1992) or the best-fit FP coordinates de- 
termined by observations in various passbands. The primary 
advantage of choosing the virial plane coordinate system is 
the easily determined tilt, which provides a gauge of possi- 
ble variations in the central total-to-stellar mass Mtotai/-^* or 
kinematic nonhomology of remnants as a function of galaxy 
mass. 

For our purposes, we choose not to use stellar population 
synthesis to compare with determinations of the FP in optical 
passbands. Observations indicate that trends in the M/L ratio 
owing to stellar population effects (e.g. age or metallicity) as 
a function of galaxy mass or luminosity will contribute sig- 
nificantly to the FP tilt, especi ally at short-wa velengths (for a 
recent result on this issue, see Cap pellari et aL.,20Q5.) . Obser- 
vations have also determined that elliptical galaxies are typi- 
cally old and their stellar populations redden passively with 
time (e.g. Bender etal. 1996; van Dokkum & Franx 1996). 
These constraints imply that to properly recover the short- 
wavelength photometric FP, precise information on the stellar 
age, metallicity, and formation-redshift distribution of ellip- 
tical galaxies as a function of stellar mass at z = must be 
obtained. We note that knowing only either the formation- 
redshift (e.g. the redshift of the last major gas-rich merger) 
or mean stellar age of elliptical galaxies may not be suffi- 
cient to determine the photometric FP of the entire ellipti- 
cal galaxy population. Recent surveys indicate that ellipti- 
cals typically undergo a major dissipationless merger at red- 
shifts z < 1 ( Bell et al. 2005; van Dokkum 2005), and, in prin- 
ciple, such mergers may induce galaxy mass-dependent tilt 
from structural effects that are disjoint from M/L effects from 
stellar populations that characteristically predate those events. 
Clearly, applying stellar population synthesis models to el- 
liptical galaxies produced in individual galaxy merger sim- 
ulations to produce simulated photometric FP scalings with- 
out attempting to correct for the cosmologically-determined 
properties of the real elliptical population is likely too naive. 
Furthermore, comparing directly stellar-mass FP scalings de- 
termined from simulations with short-wavelength (e.g. B- 
band or Sloan g-band) photometric FP scalings should be per- 
formed with extreme caution as the short-wavelength FP may 



have additional sources of tilt not present in the stellar-mass 
FP 

With these concerns in mind, our stellar-mass FP results 
wi ll be compared wi th the near-infr ared (IR) FP dete rmined 
bv lPahreetaP il998b). While, as Pahre et al. ( 1998b) note, 
tilt owing to stellar population effects may still be present 
and is not tightly constrained, structural or dynamical non- 
homologies can contribute significantly to the FP scaling in 
the near-IR. The A'-band magnitudes of stellar populations of 
a given age are less influenced by metallicity than their opti- 
cal passband luminosities (e.g. Bruzual & Chariot 2003), and 
typical /iT-band mass-to-light ratios are within ss 35% of unity 
for stellar populations formed at redshifts 0.75 < z < 5.0. 
These properties make comparisons between the stellar-mass 
and near-IR FP scalings somewhat more sensible than com- 
parisons with shorter-wavelength photometric FP scalings, 
though not ideal. 

In principle, a more straight-forward comparison would be 
to use high-resolution cosmological simulations of galaxy for- 
mation to probe the fundamental plane with simultaneously 
accounting for metallicity and stellar age effects. Previous at- 
tempts to compare results of cosmological simulations for the 
scaling la ws of elliptic al galaxies with observations have been 
made (e.g. lSaiz et al.1 '2004). but with spatial resolution some 
« 20 times worse than the isolated merger simulations pre- 
sented here. Such spatial resolution is larger than the effective 
radii of even moderately-sized (Mr > 21) early-type galaxies. 
However, the metallicities and ages of the stellar populations 
of galaxies produced in cosmological simulations can be es- 
timated throughout their formation, and the results of stellar 
population synthesis modeling would therefore be easy to in- 
terpret. 

We mention here that semi-analytic techniques could 
be combined with the results of high-resolution merger 
simulations in an attempt to account for the redshift- 
dependent formation of _the e llipti cal galaxy population. 
For example, [Robertson et al.l ll2005ah used the results of 
iHopkins et alJ j2005ah . who inferred a redshift-dependent 
black hole mass function from the quasar luminosity func- 
tion, to determine the influence of the redshift-dependent 
formation times of elliptical galaxies on the Mrh-ct rela- 
tion. Using the resu l ts from Robertson et al. (2005a) and 
IHopkins et alJ IHopkins et al. (,2005b) combined the 

redshift-dependent properties of galaxy remnants, the Mbh-ct 
relation, and the redshift-dependent black hole mass function 
to model the evolution of the red-galaxy luminosity function 
and color-magni tude relation. A combination of the results 
from IHopkins et al. ( 2005b,) and the stellar-mass FP relation 
presented here could be used to account for the effects of 
color-magnitude evolution on the photometric FP, but would 
likely involve other assumptions beyond those employed here 
and we defer such analysis for future work. 

4. RESULTS 

Below, we present the FP and R^-Mi, relations for the dis- 
sipationless, dissipational, full-model, and spheroidal merger 
simulations. For each FP and R^-M^, relation, we list the 
best fit scalings in Table|3]with \-a errors on their power-law 
slopes (i.e. A for the FP and for the R^-Mi,). 

4.1. Dissipationless Disk Simulations 

The merging of the dissipationless disks described in ^3.11 
produces a FP relation similar to the virial scaling. Figure 
n] shows the location of dissipationless merger remnants in 
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Fig. 1. — Fundamental Plane (FP) relation produced by the merg- 
ing of dissipationless disk galaxy models appropriate for redshifts 
z = (black), z = 2 (red), z = 3 (blue), and z = 6 (green) on nearly 
radial, parabolic orbits. All models include dark matter halos. The 
dissipationless merging of pure disk models (solid triangles) and disk 
models with bulges (solid circles) produce similar FP relations nearly 
parallel to the plane defined by the virial relation. Increasing the an- 
gular momentum of the orbit by lengthening the pericentric passage 
distance of the orbit produces an offset in the FP by increasing the 
effective radius of the remnants (open circles), but the systems still 
obtain a FP scaling similar to the virial plane. Select higher reso- 
lution runs closely follow the FP delineated by their lower resolu- 
tion counter-parts (open diamonds). For comparison, the best least- 
squares fit to the FP of pure disk merger remnants is plotted (dotted 
line). Also shown is the mean deviation induced by line-of-sight 
variations in projected quantities for a given remnant (detached error 
bars). 



the virial coordinate system, and plots remnants from pro- 
genitors appropriate for various redshifts. The best-fit FP 
scalings produced by pure disk systems (solid triangles) are 
a = 2.00, f3 = 1.01, with almost no tilt relative to the virial 
plane (A = 1.00). Individually, progenitors at each of the 
four simulated redshifts produce merger remnants also closely 
aligned with the virial plane (A = 0.97-1.03). Including stel- 
lar bulges in the progenitors (solid circles) produces a similar 
FP scaling {a = 1.95, /3 = 0.98; A = 0.97). A more substantial 
change in the FP is realized by increasing the pericentric pas- 
sage distance of the encounters from 2rd to 0.4/? vii (open cir- 
cles), which correspondingly increases the total angular mo- 
menta of the merging systems. Remnants produced in the dis- 
sipationless wide orbit mergers typically have larger effective 
radii, inducing an offset of Alog7?e ~ -0.4 in the FP. How- 
ever, the wide orbit FP scalings (a = 1 .97, (3 = 1 .04; A = 0.97) 
are still very similar to both the nearly radial orbit FP and 
virial scalings. Increasing the number of dark matter particles 
per halo to 180,000 in the progenitors (open diamonds) has 
very little effect on the resulting FP, suggesting our results are 
not strongly influenced by our numerical resolution. 
The R^-Mi, relation produced by merging dissipationless 



10 



Stellar Mass 

A No Bulge, z=0 | 
A No Bulge, z^2 { 
A No Bulge, z=3 < 
A No Bulge, z=S 1 
_ O Wide Orbit, z=0 < 
O Wide Orbit, z=2 
O Wide Orbit, z=3 
O Wide Orbit, z-6 
O HIgti Reeolution, z-0 



average los variation 

t Bulge, z=0 

t Bulge, z=2 

t Bulge. z=3 

t Bulge, z=6 



o 

O £y 



o 



o 

if 



O f;« 

Sji A 



Best Fit, No Bulge 
Shen et al. (2003) 



10" 



10" 



M 



10" 
[M ] 



10^' 



Fig. 2. — Effective radius ~ stellar mass M* relation produced 
by the merging of dissipationless disk galaxy models appropriate for 
redshifts z = (black), z = 2 (red), z = 3 (blue), and z = 6 (green) 
on nearly radial, parabolic orbits. All models include dark matter 
halos. The dissipationless merging of pure disk models (solid tri- 
angles) and disk models with bulges (solid circles) produce R^—Mi, 
relations shallower than that measured for massive galaxies in the 
Sloan Digital Sky Survey (Shenetal. 2003). For comparison, the 
best least-squares fit to the R^ -M* relation of pure disk merger rem- 
nants is plotted (dotted line). Also shown is the mean deviation in- 
duced by line-of-sight variations in projected quantities for a given 
remnant (detached error bars). 



disk progenitors, shown in Figure |2] reflects the trends ap- 
parent in the FP those mergers generate and additional fea- 
tures owing to redshift-dependent progenitor properties. The 
merging of pure disk galaxies (solid triangles) appropriate 
for various redshifts generates remnants with a shallower 
mean Ri^-Mi, relation {pL « 0.45, solid line) than that mea- 
sured in late-type galaxies in the Sloan Digital Sky Survey 
JShen et alJl2003L « 0.56). The remnants also systemati- 
cally decrease in effective radius with redshift, reflecting the 
decrease in progenitor disk scalelength and pericentric pas- 
sage distance. Wider orbit mergers (open circles) with larger 
angular momenta produce remnants with larger effective radii 
while galaxies containing bulges (solid circles), and therefore 
less specific angular momenta, produce remnants with smaller 
effective radii. These results are consistent with expectations 
for dissipationless systems where the energy and angular mo- 
menta of the model system are manifestly conserved. As a 
final note for the dissipationless runs, increasing the numeri- 
cal resolution of the dark matter by 3 x (open diamonds) gives 
consistent results, suggesting that artificial heating of the stel- 
lar component by discreteness effects in the dark matter halo 
is not strongly influencing the structure of the remnants. 

4.2. Dissipational Disk Simulations 

The merging of gas-rich disk galaxy progenitors including 
dissipation results in FP and R^-Mi, relations that differ sub- 
stantially from the analogous relations produced by the merg- 
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Fig. 3. — Fundamental Plane (FP) relation produced by the merg- 
ing of gas-rich disk galaxies with dark matter halos, star formation 
and supernova feedback. Shown are remnants produced by mergers 
appropriate for redshifts z = (black circles), z = 2 (red diamonds), 
z = 3 (blue triangles), and z = 6 (green squares) with nearly radial, 
parabolic orbits. The disslpational merging of pure disk models pro- 
duces a FP nearly parallel to the observed infrared FP jPahre et'SI 
^98b) and is almost independent of the redshift scalings of the pro- 
genitor systems. For comparison, the best least-squares fit to the 
FP delineated by the remnants is plotted (solid line). Also shown is 
the mean deviation induced by line-of-sight variations in projected 
quantities for a given remnant (detached error bars). 
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Fig. 4. — Effective radius ~ stellar mass M* relation produced by 
the merging of gas-rich disk galaxies with dark matter halos, star for- 
mation and supernova feedback. Shown are remnants produced by 
mergers appropriate for redshifts z = (black circles), z = 2 (red dia- 
monds), z = 3 (blue triangles), and z = 6 (green squares) with nearly 
radial, parabolic orbits. The disslpational merging of pure disk mod- 
els produces a -M* relation roughly parallel to that measured for 
massive galaxies in the Sloan Digital Sky Survey I Shen et al. 2003). 
For comparison, the best least-squares fit to the R^-Mi, relations 
delineated by the disslpational simulations (solid line) and dissipa- 
tionless simulations (dotted line) are plotted. Also shown is the mean 
deviation induced by line-of-sight variations in projected quantities 
for a given remnant (detached error bars). 



ing of dissipationless disks. Figure |3l shows the FP relation 
generated by the disslpational model disk galaxies appropri- 
ately scaled for various redshifts. The disslpational model 
FP displays a tilt relative to the virial plane (A = 0.8), with 
a scaling (a = 1.58, 13 = 0.80) similar to the near-IR photo- 
metric FP (Pahre et al. 1998b, a = 1.53, f3 = 0.79). The dissl- 
pational model FP includes remnants produced from progen- 
itors that vary by a factor 2x in gas fraction (/gas = 0.4,0.8), 
including either strongly (^eos = 1 -0) or weakly (^eos = 0.25) 
pressurized ISM equations-of-state. For gas-rich systems 
(/gas > 0.4), the gas fraction and ISM physics of the progeni- 
tors have only a slight effect on the FP. The redshift-dependent 
properties of progenitor systems have little effect on the sim- 
ulated FP plane, even as the structural properties of the disk 
galaxies vary substantially with redshift. Combined, varying 
the gas fraction above /gas > 0.4, changing the ISM pressur- 
ization dramatically, and scaUng the progenitor systems for 
redshifts z = 0-6 produce only a small amount of scatter in 

the disslpational model FP ({A^Y^^ = 0.065). 

The Re -Mi, relation of the remnants, shown in Figure |2 
is more strongly influenced by the range of progenitor prop- 
erties. The smaller progenitor galaxies appropriate for higher 
redshifts produce smaller remnants, while the less-pressurized 
ISM models also decrease the effective radii of the rem- 
nants. The disslpational model R^-Mi, relation is significantly 
steeper (solid line, /i = 0.57) than the relation produced by 



the dissipationless merging of disk galaxies (dotted line, /i = 
0.45), and compares well with the relation measured for late- 
type galaxies in the SDSS (Shenetal. 2003). The disslpa- 
tional simulations produce an R^-Mi, relation that has a lower 
/?e normalization than the Shen et al. ( 2003) relation, but sub- 
sequent re-merging and a cosmologically-representative dis- 
tribution of orbits will likely decrease the discrepancy by in- 
creasing Rg at a given stellar ma ss. Since these t he dissipa- 
tionless remnants lie above the 'Shen et al.' ("2003 ) relation, 
re-merging and higher angular momentum orbits will only 
further increase this discrepancy with real elliptical galaxies. 
While a proper accounting for the distribution of formation 
redshifts for elliptical galaxies would alter the R^-Mi, rela- 
tion by weighting the remnants unequally, the mergers of disk 
galaxies where dissipation is important will nevertheless pro- 
duce a /?e -Af* relation that is steeper than that produced by 
dissipationless merging. 

4.3. Full-Model Disk Simulations 

Introducing the effects of black hole feedback through the 
full-model simulations produces a set of fundamental scal- 
ing relations similar to those in disslpational mergers with- 
out black holes. The fundamental plane produced by the 
full-model simulations yields nearly the same FP scalings 
(a = 1 .55, (3 = 0.82) and tilt (A = 0.79) as the disslpational sim- 
ulations (see ^4.21 . and is similar to the observed near-IR seal- 
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Fig. 5. — Fundamental Plane (FP) relation produced by the merg- 
ing of gas-rich disk galaxies with dark matter halos, star formation, 
supernova feedback, and a prescription for feedback from accreting 
supermassive black holes. Shown are remnants produced by mergers 
appropriate for redshifts z = (black circles), z = 2 (red diamonds), 
z = 3 (blue triangles), and z = 6 (green squares) with nearly radial, 
parabolic orbits. The merging of pure disk galaxies using our full 
physical model produces a FP nearly parallel to the observed in- 
frared FP I Pahre et al. 1998b) and nearly coincident with the FP pro- 
duced by similar simulations without black holes. The FP relation 
is roughly independent of the redshift scalings of the progenitor sys- 
tems and the location of remnants within FP is fairly insensitive to a 
large variety of disk orientations (purple pentagons) and orbital con- 
figurations (orange hexagons), as changes in the effective radius are 
compensated by changes in the velocity dispersion and surface mass 
density. For comparison, the best least-squares fit to the FP delin- 
eated by the remnants is plotted (solid line). Also shown is the mean 
deviation induced by line-of-sight variations in projected quantities 
for a given remnant (detached error bars). 

ings dPahr e et alJl998bl) . These simulated remnants exhibit a 

1 /2 

scatter about their mean FP of (A^) ' = 0.062, comparable 
to bot h the obser ved scatter in the FP at optical wavelengths 
(e.g. Bernard! et al 2003a ) and that produced by the dissipa- 
tlonal simulations. Black hole feedback causes the full-model 
remnants to be slightly larger than the dissipatlonal model 
remnants as feedback-driven winds remove gaseous material 
from the innermost regions of the remnants that would other- 
wise contribute to the central stellar content. Figure [S] shows 
the full-model remnants produced by mergers appropriate for 
various redshifts on nearly radial, parabolic orbits. The resul- 
tant FP relation is roughly independent of the redshift scalings 
of the progenitor systems and the location of remnants within 
the FP is fairly insensitive to the orientations of the disks (pur- 
ple pentagons, see Table |2j and orbital configuration (orange 
hexagons, see Table |2ji. As these simulations indicate, vary- 
ing the total angular momentum by changing either the ori- 
entation of the disks (spin) or the pericentric passage distance 
(orbital) does not strongly influence the FP of the remnants. 
Typically the angular momentum of the orbit only influences 
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Fig. 6. — Effective radius R^^ ~ stellar mass relation produced 
by the merging of gas-rich disk galaxies with dark matter halos, star 
formation, supernova feedback, and a prescription for feedback from 
accreting supermassive black holes. Shown are remnants produced 
by mergers appropriate for redshifts z = (black circles), z = 2 (red 
diamonds), z = 3 (blue triangles), and z = 6 (green squares) with 
nearly radial, parabolic orbits. Also plotted are remnants produced 
by varying the system angular momentum through the initial disk 
orientation (purple pentagons) and pericentric passage distances (or- 
ange hexagons) for a single pair of progenitor models, which pro- 
duces a spread in the remnant effective radius. The merging of pure 
disk galaxies using our full physical model produces a R^—Mi, re- 
lation roughly parallel to th at measured for massive galaxies in the 
Sloan Digital Sky Survey I Shen et alj|2003ft . though with an offset. 
The resultant R^ -M* relation shifts downward with the redshift of 
the progenitor systems as the remnants decrease with size. For com- 
parison, the best least-squares fit to the R^-M^ relation delineated 
by these remnants (solid line), as well as the dissipatlonless disk 
model Rc-Mt relation (dotted line) and the Shen et al. (2003) re- 
lation (dashed line) are plotted. Also shown is the mean deviation 
induced by line-of-sight variations in projected quantities for a given 
remnant (detached error bars). 



the effective radii of the remnants, with the velocity dispersion 
and surface mass density changing in a compensating fashion 
to maintain the FP scalings. In addition, varying the concen- 
tration of the dark matter halo by a factor of ^ 2.5 does little 
to change the FP scalings (cyan triangles). 

The influence of angular momentum on the remnant proper- 
ties is plainly visible in the Rg-M^, relation for the full-model 
simulations (Figure|S). As is evident in the FP scalings, sim- 
ulated mergers appropriate for redshifts z = - 6 produce re- 
sults very similar to the dissipatlonal runs with a compara- 
ble R^-Mi, scaling (/i = 0.57) and slightly larger remnants. 
However, varying the angular momentum of the merging sys- 
tem through changing the progenitor disk orientations (purple 
pentagons) and pericentric passage distances for several ori- 
entations (orange hexagons) can as much as double the av- 
erage effective radius of the remnant produced by the same 
progenitors. The dark matter halo concentration also has a 
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Fig. 7. — Fundamental Plane (FP) relation produced by the merg- 
ing of spheroidal galaxy models. Using |Hernauist ( I99G) stellar 
spheroids models with dark matter halos as initial conditions (blue 
triangles), the remnants produced by dissipationless spheroid merg- 
ers appropriate for redshift z = 0, with nearly radial, parabolic orbits 
(black circles) and wide, elliptical orbits with circularity e = 0.5 (red 
diamonds) are calculated. The remnants of spheroidal mergers pro- 
duce a FP relation similar to their progenitor systems, roughly inde- 
pendent of the orbital energy or angular momentum. Also shown are 
select spheroidal remnants from gas-rich disk galaxy merger simu- 
lations that include star formation, supernova feedback, and a pre- 
scription for feedback from accreting supermassive black holes (or- 
ange squares). The disk galaxy merger remnants occupy a FP re- 
lation similar to that observed in infrared observations (solid line). 
The re-merging of these disk galaxy remnants on nearly radial, 
parabolic orbits (green circles), further demonstrates that the merg- 
ing of spheroidal remnants produ ces a FP similar to t hat occupied 
by the progenitor systems (e.g. Caoelato et al. 1995; ! Dantas et alj 
1^03; NiDoti et al. 2003; Bovlan-Kolchin et al. 2005). For compar- 
ison, the best least-squares fit to the FP delineated by the dissipa- 
tionless disk merger remnants (dotted line) and spheroidal merger 
remnants (dash-dotted line) is plotted. Also shown is the mean devi- 
ation induced by line-of-sight variations in projected quantities for a 
given remnant (detached error bars). 
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Fig. 8. — Effective radius - stellar mass M* relations pro- 
duced by t he merging of spheroidal galaxy models. The mergers of 
iHemguisil fT990i') stellar spheroids initial conditions with dark mat- 
ter halos (blue triangles), initialized to the Shen et al. (2003) rela- 
tion (dashed line), with nearly radial, parabolic orbits (black circles) 
and wide, elliptical orbits with circularity e = 0.5 (red diamonds) are 
simulated. The remnants of spheroidal mergers produce a R^-Mi, 
relation slightly shallower than the observed relation, roughly inde- 
pendent of the orbital energy or angular momentum. Also shown 
are spheroidal remnants from gas-rich disk galaxy merger simula- 
tions that include star formation, supernova feedback, and a prescrip- 
tion for feedback from accreting supermassive black holes (orange 
squares). The disk galaxy merger remnants occupy a -M* rela- 
tion with the same scaling as observed in SDSS (dashed line). The 
re-merging of these disk galaxy merger remnants on nearly radial, 
parabolic orbits (green circles) produces a similar R^-Mi, relation 
as the dissipationless merging of spheroidal galaxies. For compari- 
son, the best least-squares fit to the R^ -M* relation delineated by the 
dissipationl ess sp heroidal merger remnants (dashed-dotted line), the 
IShen et alj <2003t) relation, and the relation produced by simulations 
using the full physical model (solid line) are plotted. Also shown is 
the mean deviation induced by line-of-sight variations in projected 
quantities for a given remnant (detached error bars). 



noticeable effect on the galaxy size (cyan triangles). We reit- 
erate that these changes in the effective radius only induce a 
spread parallel to the FP (see above). Moreover, as is the case 
for the dissipational simulations without black hole growth, 
the R^-Mi, relation for the full model remnants would only 
benefit from subsequent re-merging or more cosmologically 
representative orbits by increasing R^. 

4.4. Spheroid Simulations 

As demonstrated in ^4.H - R31 the merging of dissipation- 
less disk galaxies produces a FP with almost no tilt relative to 
the virial plane and a R^-Mi, relation that is shallower than 
observed, while gas-rich disk galaxies yield FP relations with 
substantial tilt and reasonable R^-Mi, scalings. A natural ex- 
tension of these calculations is to test whether the subsequent 
merging of the remnants or model galaxies similar in struc- 



ture to the remnants preserves the FP and R^—M^, relations 
produced by the first generation of merg ers. The spheroidal 
galaxy model progenitors described in ^3.4 | that are initial- 
ized to obey the observed Shen et al. ( 2003) Rg -M* relation 
delineate a tight FP relation (blue triangles, Figure^. Merg- 
ing these spheroidal galaxy models on either nearly-radial 
parabolic orbits (red diamonds), wide slightly-bound orbits 
(black circles), or in an exactly head-on encounter (purple 
pentagon) produces a similar tight FP relation, only slightly 
tilted relative to the progenitor FP scalings. As described pre- 
viously, altering the angular momentum of the orbits has the 
effect of increasing the effective radii of the remnants but does 
not strongly affect the FP relation as the velocity dispersion 
and surface mass density compensate to roughly maintain the 
progenitor FP scalings. 

A similar result is obtained by re-merging remnants from 
the full-model simulations. The FP delineated by the full- 
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model remnants (Figure solid line) is roughly recovered 
by re-merging the disk galaxy remnants, which demonstrates 
that the approximate preservation of the FP relation holds for 
relations instilled by hand (as for spheroidal model galax- 
ies) and those generated self-consistently through the merg- 
ing process (as for the full-model merger remnants). The 
re-merger remnants are consistently larger in effective radii 
when compared with their progenitors, also in agreement with 
the spheroidal galaxy mode l mergers and previous work (e.g. 
iBovlan-Ko lchin et alJ2005h . 

The behavior of the R^-Mi, relation under merging demon- 
strates the increase of R^ of spheroidal and re-merger rem- 
nants relative to their progenitor p roperties . As noted previ- 
ouslv jBarnesll!9 92: Hernauist et a l.l|1993t IC ole et al. 2003 
Nipoti et all l2003t iShen et al. 2003; 'Bovlan-Kolchin et 3 
the effective radius of the spheroidal merger remnants 
approximately doubles in an equal-mass merger and will in- 
troduce scatter parallel to the R^-ax\s. However, as is apparent 
from Figure |8] the R^-Mi, relation only moderately changes 
as a mass-sequence (e.g. ^ = 0.56 for spheroidal progenitors, 
^ = 0.52 for spheroidal merger remnants). Re-mergers of disk 
galaxy remnants produce a similar change in the R^-M^, re- 
lation. A dissipationless merger between spheroids then is 
expected to introduce scatter into the R^ - relation but 
does not completely destroy the /?e - correlation of the 
progenitors. The presence of a well-defined R^-Mi, corre- 
lation at z = is therefore not a strong constraint on the dis- 
sipationless merging of spheroids (c.f. IShen et al.ll200l . al- 
though substantial additional merging beyond the recently in- 
ferred 0.5-1 dissipationless mergers at redshifts z < 1 typical 
for spheroidal galaxies (Bell et al. 2005; van Dokkum 2005} 
could possibly induce noticea ble scatter into th e R^-Mi, rela- 
tion beyond that reported by S hen et alJ (120031) . 

We note briefly here that the maintenance of a R^-Mi, re- 
lation with ji K, 0.5 also will keep the Mbh-c relation scaling 
of Mbh oc a'^, while possibly inducing some additional scat- 
ter owing to a shift in the normalization toward lower velocity 
dispersions at a given black hole mass. The velocity disper- 
sion of a spheroidal galaxy is roughly 



The maintenance of the Re—Mi, relation with /i sa 0.52 after 
re-merging and the assumption that the stellar mass and 
black hole mass Mbh both double then implies 

Mbh (X (74-2^ (18) 

which is approximately the observed Mbh-c scaling at z = 0. 
Correspondingly, the Mbh-c relation of the re-merger rem- 
nants that include black holes (see Figure |5} does not dif- 
fer strongly from the Mbh-c relation delineated by their 
progenitor systems. While a does not increase dramat- 
ically enough to perfectly preserve the observed relation 
dNipoti et al.ll2003h . the measured Re -Mi, relation of the re- 
merger remnants provides a preliminary explanation for the 
rough maintenance of the Mbh-c relation through infrequent 
dissipationless mergers. We also note that if dissipation- 
less major merging happens more frequently for more mas- 
sive systems, then the Mbh-o" relation may also steepen rel- 
ative to the Tremaine et al. (2002) best-fit relation at higher 
masses. A signature of such dynamical reprocessing may ap- 

S iear as a mass-depende nt scatter in the local Mbh-<7 relation 
Robertson et alJ2005bl) . 
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Fig. 9. — Black hole mass Mbh ~ stellar velocity dispersion relation 
a relation measurements for galaxies produced by the re-merging 
of disk galaxy merger remnants that obey the local Mbh-o" rela- 
tion (e.g., 2'remaine et al. 2002, solid line). While the number of 
re-mergers examined is limited, a single generation of dissipation- 
less mergers after the initial formative gas-rich mergers that generate 
the Mbh-o" relation is not expected to strongly alter the observed 
Mbh-o" relation but may be a source of scatter. 



5. THE IMPACT OF PROGENITOR GAS FRACTION ON MERGER 
REMNANT PROPERTIES 

The results presented in ^^demonstrate that the merger of 
gas-rich systems produces a substantially different fundamen- 
tal plane scaling than the virial scalings obtained by the merg- 
ing of dissipationless disk galaxies. However, the origin of the 
FP tilt measured for the simulations that include the effects of 
dissipational physics has not yet been illuminated. To better 
demonstrate that dissipational physics does indeed produce a 
FP tilt for sufficiently gas-rich systems and to demonstrate the 
manner by which gas dissipation gives rise to FP tilt, we now 
present the results of a sequence of dissipational model merg- 
ers in which the gas fraction /gas of the progenitor disks has 
been varied systematically from /gas = 0.01 t o /ga s = 0.4. The 
details of these simulations are described in ^3.2l and listed in 
TablelU 

5.1. Fundamental Scaling Relations 

The Fundamental Plane relation produced by the merging 
of disk galaxies models with dissipation as a function of pro- 
genitor gas content is shown in Figure^] The remnants pro- 
duced by progenitor systems with various gas fractions, /gas 
are plotted. The impact of dissipation on the FP tilt of rem- 
nants exhibits a behavior where for remnants of progenitors 
with /gas < 0.2 little FP tilt is induced by gas dissipation while 
for remnants of progenitors with /gas > 0.2 substantial FP tilt 
is induced (see Table^for the best-fit scalings). The nearly- 
dissipationless simulations (/gas < 0.1) produce FP scalings 
similar to the virial scalings (A = 0.92-0.97), while the gas- 
rich mergers (/gas = 0.4) produce a substantial tilt (A = 0.83). 
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The Re -M^ relation of the remnants also displays a behav- 
ior where above a critical gas fraction /gas « 0.3 dissipation 
is important for determining the remnant properties. Below 



0.2 t he Rg-M^ relation is shallower than the observed 

IShen e t aL (2003) relation and is similar to the Rg-M^, re- 
lation produced by dissipationless merging of purely stellar 
disks. As the gas fraction of the progenitor systems increases 
beyond /gas = 0.2, the dissipational effects of gas cooling and 
star formation produce smaller remnant effective radii and a 
steeper R^-Mi, relation (see Figure [TTl. The decrease in ef- 
fective radius is most dramatic in lower-mass systems, where 
Rg can decrease by as much as 50%. The tilt induced in the FP 
relation with progenitor gas fraction by the increasing impor- 
tance of gas dissipation appears closely related to this mass- 
dependent decrease in effective radius, as the velocity disper- 
sion and surface mass density adjust with the change in effec- 
tive radius to produce FP tilt as the dissipative effects increase. 
We note here that the equation of state of the ISM also influ- 
ences the effective radii, with the less-pressurized ISM model 
leading to smaller remnants as dissipation becomes more ef- 
ficient with decreasing ISM pressurization. 

5.2. Phase-Space Density 

While the behavior of the fundamental scaling relations of 
the merger remnants as a function of progenitor gas fractions 
demonstrates that dissipation can induce FP tilt by steepening 
the Rg-M^, relation, these structural parameter correlations 
reflect the composition of more fundamental properties of the 
remnants. One such property is the fine-grained phase-space 
density f(x, v) at the spatial coordinate x and velocity coordi- 
nate V. The quantity f(x, v)(ixdv describes the probability of a 
stellar mass element dM^, being located within the differential 
phase-space volume djcdv centered at coordinates (x,v). For 
a purely collisionless system, f(x,v) satisfies the Boltzmann 
equation 

9f ...... „^df 



^+v«V/-V«'»Tr- = 0. 

dt d\ 



(19) 



Also known as the Vlasov formula. Equation ( I19> describes 
the detailed conservation of the phase-space density as an in- 
compressible flow under the influence of a smoothly varying 
potential $(x,f) and represents the 1 -particle limit of Liou- 
ville's theorem. The underlying phase-space structure of a 
dissipationless system is then a conserved, fundamental quan- 
tity. 

Although the fine-grained phase-space density of stars is a 
basic physical property of a galaxy, the finite numerical reso- 
lution of simulations limits the measurable phase-space den- 
sity to the quantity /(x,v)AxAv that represents the probabil- 
ity of a macroscopic stellar mass element AM* being located 
within the macroscopic phase-space volume AxAv centered 
at (x, v). The quantity / is commo nly referred to as the coarse- 
grained distribution function (e.g. lBinnev & Tremainal987h . 
While /(x,v) can be measured from the simulations by grid- 
ding the stellar particles into AxAv volumes, a well-sampled 
f(x,v) requires a cru de AjcAv binning. I nstead, one may 
adopt the approach of Hernaui st et alJjl993^ and measure the 
coarse-grained phase-space density f{E) averaged over hy- 
persurfaces with energy E under the assumptions of spheric- 
ity and isotropy. With these approximations, the fine grained 
phase-space density / is a function of energy and given a 
model of a galaxy that provides the density of states g{E), the 
function describing the phase-space volume of energy hyper- 
surfaces, one may relate the mass-energy distribution of stars 
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Fig. 10. — Fundamental Plane (FP) relation produced by the merg- 
ing of disk galaxies with dark matter halos, star formation and super- 
nova feedback as a function of progenitor gas content. The galaxy 
models are appropriate for z = and are merged on nearly radial, 
parabolic orbits. Shown are remnants produced by progenitor sys- 
tems with gas fraction /gas = 0.01 (red squares and line), /as = 0.025 
(yellow triangles and line), /as = 0.04 (green diamonds and line), 
/as = 0.1 (blue triangles and line), /as = 0.2 (purple hexagons and 
line), and /as = 0.4 (black circles and line). As the gas fraction of 
the progenitor systems increases past /as = 20%, the dissipational ef- 
fects of gas cooling and star formation induce a substantial tilt in the 
resulting FP relation. Nearly dissipationless progenitor (/as < 20%) 
produce FP relations roughly parallel to the plane delineated by the 
virial relation, while for systems where gas dissipation is important 
the res ultant FP is similar to the observed infrared FP ( Pahre et aO 
119989). The increasing tilt is generated through a decrease in the 
remnant effective radius with increasing gas content and the corre- 
sponding change in the central stellar mass fraction (see also Figures 

QII-IIU. 



dM*/ dE to the fine-grained phase-space density /(£) as 



dM* 



(20) 



A coarse grained distribution function f(E) can be measured 
from the simulations through a histogram of stellar particle 
energies AM^,/ AE as 

1 AM* 

/(£)=— (21) 
giE) AE 

As iHernauist et alJ Jl993h demonstrate, the coarse-grained 
/(£■) provides a reUable estimate of the fine-grained /(£) for 
numerical models of isotropic, spherical galaxies. An even 
more precisely measurable quantity is the cumulative distri- 
bution function 

Jeo dE 

where Eq is the central binding energy. The cumulative dis- 
tribution function s(f) provides the fraction of stellar mass 
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Fig. 11. — Effective radius i?c - stellar mass M* relation produced 
by the merging of gas-rich disk galaxies with dark matter halos, star 
formation and supernova feedback as a function of progenitor gas 
content. The galaxy models are appropriate for z = and are merged 
on nearly radial, parabolic orbits. Shown are remnants produced 
by progenitor systems with gas fraction /gas = 0.01 (red squares and 
line), /gas = 0.025 (yellow triangles and line), /gas = 0.04 (green dia- 
monds and line), /gas = 0. 1 (blue triangles and line), /gas = 0.2 (purple 
hexagons and line), and /gas = 0.4 (black circles and line). As the gas 
fraction of the progenitor systems increases, the dissipational effects 
of gas cooling and star formation produce smaller remnant effective 
radii and a corresponding change in the central stellar mass fraction 
(see also Fieures fTolll4l and ll5> . These dissipational effects gener- 
ate a substantial tilt in the Fundamental Plane above a gas fraction of 
roughly /gas > 30%. 
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Fig. 12. — Cumulative coarse-grained distribution function i(/) of 
a remnant from the equal-mass merger of progenitors with Kir = 115 
km s"' as a function of progenitor gas fraction. Shown is the frac- 
tion of remnant stellar mass with phase-space densities greater than 
/ for mergers with /gas = 0.01 (red line), /gas = 0.025 (orange line), 
/gas = 0.05 (green hne), /gas = 0.1 (blue hne), /gas = 0.2 (purple line), 
and /gas = 0.4 (black line). For gas fractions /gas < 0.2, the phase- 
space density distributions of the remnants do not change apprecia- 
bly; the binding energies of the stars do not increase greatly relative 
to those in a collisionless merger as dissipation is not efficient at al- 
tering the structure of the galaxy. For gas fraction /gas > 0.2, the 
phase-space density distribution of the remnant shifts toward higher 
values and more strongly-bound energy hypersurfaces. For low mass 
ellipticals the phase-space density of remnants is then increased in 
gas-rich /gas > 30% mergers where dissipation is important for the 
remnant properties. 



at phase-space densities above /, which is equivalent to the 
fraction of stars more tightly bound than energy E(f) when 
averaged over energy hypersurfaces. For a binned energy his- 
togram AM^/A/i measured from the simulations, a cumula- 
tive coarse grained distribution function s{f) can be calculated 
as 

-^^J_g(AM^^ (23) 

where (AM^),- is the stellar mass with energy between £, and 
Ei + AE, = is the index of the energy bin with E =Eq and 
/^(y) is the index of the energy bin with energy E = E(f). The 
cumulative energy distribution s{E) can be measured directly 
from the simulations, independent of any assumed model for 
g{E). The choice of galaxy model then simply correlates / 
with an energy E. 

To connect a coarse-grained phase-space density / with 
an energy E, we adopt an isotropic, stellar |Hernauist ( 1990) 

fheroid embedded within a larger isotropic, dark matter 
5rnauist (1990) spheroid (hereafter the HH model). Ciotti 
996) analytically calculated both the distribution function 
/(£) and density of states g{E) for such a system, which 
serves as a convenient galaxy model for the remnants and pro- 



vides a conversion between binding energy and phase-space 
density for a set of remnant scale radii a* and adm and stellar 
and dark matter masses. The scale radii of the remnants can be 
related to the effective radi us an d reliably estimated from the 
particle distributions (e.g. ^3.5> . while the mass of the stellar 
and d ark matter components are tracked explicitly. However, 
as the iHernauisll (.1990) model may not perfectly reflect the 
stellar distribution, the potential well depth at the center of 
the merger remnant often exceeds the expected potential min- 
imum for the HH model. In this case we find the effective 
cii, for the HH system that produces the coiTect minimum po- 
tential and use this model to calculate f(E) and g{E). As the 
functional form of f(E) for the HH model depends only on 
the ratios M^^/Mdm and a^/flDM this technique for measuring 
s{f) is fairly insensitive to the manner in which we measure 
galaxy properties. Once violent relaxation reduces the time- 
dependence of the galaxy potential during a dissipationless 
merger, the form ofsif) will be fixed, while for gaseous merg- 
ers, the form of s{ f) will stabilize after the additional process 
of star formation becomes relatively unimportant. A compar- 
ison of the_5(/) resulting from strongly dissipative mergers 
with the s{f) of dissipationless remnants will indicate the im- 
pact of dissipation on the final galaxy phase-space structure. 
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Fig. 13. — Cumulative coarse-grained distribution function of 
a remnant from the equal-mass merger of progenitors with Vvir = 500 
km s"' as a function of progenitor gas fraction. Shown is the fraction 
s{J) of remnant stellar mass with phase-space densities greater than 
/ for mergers with /gas = 0.01 (red line), /gas = 0.025 (orange line), 
/gas = 0.05 (green line), /gas = 0.1 (blue hne), /gas = 0.2 (purple line), 
and /gas = 0.4 (black line). For gas fractions /gas < 0.2, the phase- 
space density distributions of the remnants do not change apprecia- 
bly; the binding energies of the stars do not increase greatly relative 
to those in a collisionless merger as dissipation is not efficient at al- 
tering the structure of the galaxy. For a gas fraction /gas > 0.2 the 
phase-space density distribution of the remnant increases at mod- 
erate phase-space densities relative to the dissipationless remnant, 
but not dramatically. The comparatively small shift in s(f) for mas- 
sive ellipticals formed in gas-rich mergers relative to the larger re- 
sultant shift in low mass elliptical s results in a mass-dependent trend 
inM,otai/M* (see FiguresfTil-fTsl. 



Figure[T2lshows the cumulative coarse-grained distribution 
function ,?(/) of a remnant from the equal-mass merger of 
progenitors with Vvii = 115 km s"' _as a function of progeni- 
tor gas fraction. The measured s(f) for gas fractions below 
/gas < 0.2 remains roughly constant. The binding energies of 
the stars are not strongly affected by the weak influence of dis- 
sipation on the galaxy structure for low gas fraction merger, 
but as the progenitor gas fractions increase beyond /gas > 0.2 
the phase-space density distribution of the remnant shifts to- 
ward higher values and more strongly-bound energy hyper- 
surfaces. Strongly dissipational mergers clearly lead to higher 
phase-space densities in low mass elliptical galaxy remnants. 
We note here that the shapes of the cumulative phase-space 
density as a function of mass and gas fraction remain very 
similar for either ISM pressurization we study in this work 
(?EOS = 0.25 and ^eos = 1.0), except for a slight decrease in 
the maximum phase-space density reached in the ^eos = 10 
simulations. 

For the impact of dissipation to affect the FP tilt or R^-M^ 
slope, the degree to which the phase-space density increases 
during gas-rich mergers must vary as a function of galaxy 
mass. Figure shows the cumulative coarse grained dis- 



FiG. 14. — Ratio of total mass Mtoiai to stellar mass M* within a 
spherical radius for the remnants of equal mass mergers of = 115 
km s"' disk galaxies with dark matter halos, star formation and super- 
nova feedback as a function of progenitor gas content. The galaxy 
models are appropriate for z = and are merged on nearly radial, 
parabolic orbits. Shown is Mtotai/A^* vs. R/R^ for remnants pro- 
duced by progenitor systems with gas fraction /gas = 0.01 (red line), 
/gas = 0.025 (yellow line), /gas = 0.04 (green line), /gas = 0.1 (blue 
line), /gas = 0.2 (purple line), and /gas = 0.4 (black line). For small- 
mass spheroidal remnants, as the progenitor gas content is increased 
the stellar content substantially increases in the very central regions 
(R < 0.3Rc) of the galaxy. The decrease in Mtotai/M* follows the 
corresponding decrease in i?c as star formation is concentrated in the 
high central gas overdensities induced by the merger. 



tribution function s{f) for remnants of a fiducial equal-mass 
merger between disk galaxies with Vvir = 500 km s"^ as a func- 
tion of progenitor gas fraction. While the behavior for low 
gas fractions /gas < 0.2 is similar between low mass ellipti- 
cals (Figure and massive ellipticals, the strong shift seen 
in the phase-space density distribution of low mass ellipticals 
is not measured for the more massive system. The simulations 
show that the impact of dissipative physics on the phase-space 
density of ellipticals is mass-dependent, being more important 
for low mass galaxies. 

5.3. Central Mass-to-Light Ratio 

With a measured mass-dependent change in phase-space 
density of remnants produced in gas-rich mergers, the phase- 
space density can be directly connected to changes in the FP 
tilt and R^-M^, slope induced by strongly dissipational merg- 
ing. The correspondence between phase-space density / and 
binding energy demonstrates that the regions at the center of 
the galactic potential are gaining stellar mass density. The 
gas dissipation has a less dramatic effect on the dark matter 
distribution of the galaxy, and while the gas does change the 
dark matter profile by gravitationally dragging the dark matter 
inwards, the effective radius of the dark matter decreases by 
only 250/z-' pc (0.6%) between the /ga,, = 0.1 and f„.^, = 0.4 
runs for Vvh = 115 km s"' halos. Dissipation affects the stellar 
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Fig. 15. — Ratio of total mass Mtoiai to stellar mass M* with spheri- 
cal radius for the remnants of equal mass mergers of Kir = 500 km s"' 
disk galaxies with dark matter halos, star formation and supernova 
feedback as a function of progenitor gas content. The galaxy models 
are appropriate for z = and are merged on nearly radial, parabolic 
orbits. Shown is M(otai/W* vs. R/R,: for remnants produced by pro- 
genitor systems with gas fraction /gas = 0.01 (red line), /gas = 0.025 
(yellow line), /.as = 0.04 (green line°), /ga,, = 0. 1 (blue line)° /.as = 0.2 
(purple line), and /gas = 0.4 (black line). For large-mass spheroidal 
remnants, as the progenitor gas content is increased the central stel- 
lar content begins to increase substantially after R > 0.3Rs. The 
Motai/M* value near the center is larger for massive systems than 
for less-massive spheroids (see Figure [141 . reflecting the increasing 
importance of dark matter at the center of massive spheroidal galax- 
ies. 
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Fig. 16. — Ratio of the dynamical mass estimator Mdyn = a Rc/G 
to the stellar mass M* in the full-model simulations as a function 
of stellar mass for remnants of mergers between progenitor disk 
galaxies appropriate for redshifts z = (black circles), z = 2 (red dia- 
monds), z = 3 (blue triangles), and z = 6 (green squares). Also shown 
is the measured Mdyn/M* for z = mergers with varying disk ori- 
entations (purple pentagons), orbital parameters (orange hexagons), 
and dark matter concentrations (cyan triangles). The best-fit trend 
Mdy„/M* oc '** (solid line) recovers the complete fundamental 
plane tilt measured for the full-model simulations in Figure |5| We 
also estimate the tilt contributed by variations of Mtotai/M* averaged 
within an effective radius as a function of stellar mass measur ed from 
the remnant particle distribution (shaded area, see 95.3. it . which 
contributes ~ 40- 100% of the nonhomology-related FP tilt. The 
Mtotai /A^* tilt estimates shown have been renormalized by a constant 
to intersect with the best-fit Mdyn/M* relation at the same location. 



mass much more dramatically than the dark matter mass, and 
the net result is a radius-dependent change in the Mtotai/A^* 
ratio as a function of progenitor gas fraction. 

Figure[3]shows the ratio of total mass Mtotai to stellar mass 
within a spherical radius R for the binary Vyir =115 km 
s"' mergers with various progenitor gas fractions, /gas. As the 
effects of dissipation increase, the central regions of the low- 
mass remnant become progressively more baryon-dominated. 
The decrease in the central Mtotai /A^* follows the correspond- 
ing decrease in /?e as star formation is concentrated in the high 
central gas overdensities induced by the merger. The value of 
Mtovdi/Mi, does not change appreciably at /? = /?e until the gas 
fraction reaches /gas = 0.1 and decreases dramatically as the 
gas fraction increases toward /gas = 0.4. 

The impact of dissipation on Mtotai/A^* of more massive 
systems is less important, as expected from the behavior of 
the phase-space density and 7?e of those systems. Figure [Tsl 
shows the ratio Mtotai /A^* with spherical radius R for the rem- 
nants of equal mass mergers of Vyii =500 km s" 1 systems with 
various gas fractions. For high-mass spheroidal remnants, as 
the progenitor gas content is increased, the stellar content be- 
gins to increase substantially at radii R > 0.3/?e- The cen- 
tral Mtotai /A^* remains relatively unaffected by dissipational 
physics and remains roughly constant Mtotai /M,^ ~ 1.25. The 



value of Mtotai /A^* near the center is larger for massive sys- 
tems than for less-massive spheroids (see Figure [T4t. reflect- 
ing the increasing importance of dark matter at the center of 
massive spheroidal galaxies and contributing tilt to the FP. 

5.3.1. Dynamical Mass-to-Light Ratio 

For the stellar-mass FP relation, an observationally relevant 
quantity is the estimator of the dynamical mass within an R^ 

a^R^ 

Mdyt, = k—^ (24) 

where G is the gravitational constant and A; is a constant that 
depends on the structure of the galaxy {k « 7.25 for a kinemat- 
ically isotropic Hernauist ( 1990) spheroid). The dynamical 
mass estimator allows us to characterize the FP tilt in terms of 
the observationally-accessible mass ratio 

^ocMj. (25) 
M, * 

For a stellar-mass FP of the form R^ oc a"I^^ and an R^-Mi, 
relation of the form /?e oc M* , the power-law index 7 can be 
estimated as 

7 = Ai(2/a-l-l-4/3/a) + 2/3/a-l. (26) 



Gas Dissipation & the FP 



17 



0.0 




► V =115kms' 

vir 

• V =500 km s 



2.0 



-0.2 0.0 0.2 0.4 0.6 0.8 1.0 
t-t [h-^ Gyr] 

disk 



1.2 



Fig. 17. — Mass fraction of star forming gas for Kir = 115 km s"' 
(blue triangles) and Wir = 500 km s"' (red circles) progenitor merger 
remnants in simulations with gas fraction /gas = 0.4 as a function of 
time /-?disk since the disks first interact. The gas in the massive sys- 
tem is collisionally heated during the first passage and subsequently 
cannot cool and form stars, while gas in the small mass system can 
cool sufficiently to continue to form stars throughout the remainder 
of the merger. 



For the virial scalings (a = 2, /3 = 1), the dynamical mass es- 
timator is proportional to the stellar mass and there is no FP 
tilt (7 = 0), in depe ndent of the R^-M^, relation. The results 
presented in ^.31 indicate that the full-model FP (a = 1.55, 
(3 = 0.82) is tilted relative to the virial scalings. Combined 
with the full-model -M^ relation (p. = 0.57), we expect that 
A^dyn/A^* cx M° '^. Figure [T6lshows the quantity Mdyn/M* as 
a function of stellar mass for the full-model remnants from 
mergers of progenitor disk galaxies appropriate for redshifts 
z = 0-6. While there is considerable scatter in Mdyn/M^, 
which mostly arises from the linear dependence of Mdyn on R^, 
the best fit trend in the simulations has a scaling Mdyn/M* cx 
(solid line) that recovers essentially the complete tilt in- 
ferred for the full-model FP. 

In principle, the tilt in the stellar-mass FP inferred from the 
dynamical mass estimator may include contributions from to- 
tal mass-to-stellar mass Mtotai/A^* variations, kinematic non- 
homologies that affect the velocity dispersion as a function 
of stellar mass, or other sources of nonhomology. Figure 
^]shows the estimated range of tilt contributed by Mtotai/A^* 
variations with stellar mass averaged over an effective radius 
as measured from the particle distributions of the remnants 
(shaded area), bounded by least-squares fits to the Mtotai/A/* 
trend with mass minimizing with respect to either Mtotai/A/* 
orM*. Based on this comparison, we estimate that variations 
of Mtotai/A^* with stellar mass averaged over an effective ra- 
dius contribute w 40- 100% of the nonhomology-related tilt 
in the full-model simulations. 

5.4. Dissipation vs. Galaxy Mass 
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Fig. 18. — Consumption timescale of star forming g as in Vvii' — 115 
km s"' (blue triangles) and Vyji = 500 km s"' (red circles) progeni- 
tor mergers measured for simulations with gas fraction /gas = 0.4 and 
equation of state parameter ^eos = 0.25. The consumption timescale 
is plotted as a function of time / - /disk since the disks first interact 
and reflects the ratio of Mgas(p > psfr)/(dM*/d/) during the simula- 
tion for gas above the density threshold for star formation. The cold 
gas in the massive system that is not collisi onal ly heated by the first 
pericentric passage into the halo (see Figure fTTl is always efficiently 
converted into stars. In smaller systems, the cold gas in the system 
is efficiently converted into stars only during times of interactions, 
such as pericentric passage (f-/disk = 0.1-0.2 h"' Gyr) or final coa- 
lescence (/-fdisk ~ 0.9 h"' Gyr), but is consumed on long timescales 
otherwise. 



The previous sections have demonstrated that gas dissipa- 
tion impacts the fundamental scaling properties of ellipticals 
by altering their central stellar content in a mass-dependent 
manner, being more important for small mass systems. We 
now explain how this mass-dependent impact of dissipation 
arises in our simulations. 

Figure shows the mass fraction of dense, star forming 
gas in merging disk galaxies as a function of the time t - fdisk 
since the gaseous disks first interact for a simulation with pro- 
genitor gas fraction /gas = 0.4 and equation of state parameter 
?EOS = 0.25. In massive systems, such as the Vvh- = 500 km 
s"' progenitor equal-mass merger shown here (red circles), a 
large fraction of the gas is collisionally heated during the first 
pericentric passage. This hot gas becomes x-ray luminous but 
cannot cool efficiently (Cox et al. 2004, 2005a), and mostly 
remains in the halo for the duration of the merger The effi- 
ciency of halo gas cooling decreases with halo mass, causing 
more massive systems to act increasingly like dissipationless 
systems. In small systems (we show an equal-mass merger 
of Vvir =115 km s"' progenitors, blue triangles), collisionally 
heated gas can cool efficiently and return to the disk to con- 
tinuously form stars. The continued presence of star forming 
gas in small mass systems allows for very high phase-space 
density stars to form in abundance during the height of the 
merger, decreasing the central Mtotai/A^* in the final remnant. 
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Although colhsional heating of the gas in massive systems 
limits the fraction of gas dense enough to form stars, the re- 
maining dense gas is efficiently turned into stars throughout 
the merger Figure^jshows the consumption timescale 

^consumption = -^gasCP > Psfi) 

which indicates the characteristic timescale for the mass of 
gas at densities above the star formation threshold Mg^sip > 
Psfi) to be converted into stars at the current star formation 
rate dM^/df. The massive system is always efficiently con- 
verting the available gas into stars, roughly independent of the 
progenitor separation after first approach when tidal deforma- 
tion of the disk begin vigorous star formation. We note that 
while the star formation is always efficient, the massive sys- 
tem, whose original progenitors had gas fractions of /gas = 0.4, 
becomes largely dissipationless after t - fdisk = 0.4 /i"' Gyr 
when the total gas fraction drops below /gas = 0.06 and the 
amount of cold, star forming gas is a fraction ^ 0.30 of the gas 
mass. Dissimilarly, the small mass system experiences large 
variations in its consumption timescale, with extremely effi- 
cient star formation at the first pericentric passage (f - fdisk ~ 
0.1-0.2 /i"' Gyr) and the final coalescence (f-?disk ~ 0.85 
Gyr) and inefficient star formation at larger progenitor sepa- 
rations and after the merger is completed. We comment that 
this variation in consumption timescale for small ma ss sys- 
tems d iffers importantly from the scenario described bv lBekkH 
Jl998h . who infers that tilt in the fundamental plane may arise 
from a gradation in consumption timescale or star formation 
"rapidity" with galaxy mass. While the time-averaged con- 
sumption timescale in smaller systems is longer as empha- 
sized by .BekkL (.1998:.) . we do not introduce an additional 
ad hoc mass-dependency to the gas consumption timescale. 
Differences in the gas consumption timescale as a function 
of galaxy mass in our modeling owe to the increase in disk 
gas density with galaxy mass, with low mass systems natu- 
rally having comparatively inefficient star formation relative 
to massive systems in quiescence. However, during phases 
of strong interaction, the gas density in low mass systems in- 
creases dramatically and leads to a rapid decrease in the gas 
consumption timescale. The dynamism of the consumption 
timescale in small mass systems then allows the less massive 
system to "save" its gas for strong interactions, when dissipa- 
tion is very effective at converting gas into stars in the central- 
most regions on a short timescale (e.g. Mihos & Hernquist 
(1991 . 

6. DISCUSSION 

Our simulations indicate that mergers of gas-rich disk 
galaxies produce tilt in the Fundamental Plane relationship 
relative to the scalings predicted by the virial theorem for ho- 
mologous systems. The remnants of gas-rich mergers are not 
homologous as gas dissipation is more effective at redistribut- 
ing stellar material relative to dark matter in smaller systems 
than in larger ones. Strongly dissipational mergers also appear 
necessary to steepen the R^-Mi, relation o f disk galaxy rem- 
nants to more closely follow the observed IShen et all (120031) 
relation for early-type galaxies. 

Most conservatively interpreted, our results demonstrate 
simply that if elliptical galaxies originate in gas-rich mergers 
then they will obey a tight fundamental plane relation that has 
significant tilt owing to dissipative effects including central 
Motai/A/* variations as a function of galaxy mass. However, 



independent observational and theoretical evidence suggests 
that gas-rich mergers are a necessary feature of the merger 
hypothesis as an elliptical galaxy formation scenario. 

Bender et al. (1992) provided an early characterization of 
the importance of dissipation varying with the mass of ellip- 
tical galaxies. In their "gas/stellar" (GS) continuum, massive 
systems experience less dissipation, slowly rotate, and are in- 
creasingly flattened by anisotropy. Our work describes the 
impact of the GS continuum on the central structure of el- 
liptical galaxies as a function of their stellar mass and de- 
tails the origin of the GS continuum in terms of the ineffi- 
cient cooling of coUisionally heated gas in massive systems 
and dynamic variations in the gas con sump tion timescale in 
smaller mass mergers. As we note in ^5.41 these effects dif- 
fer substantially in detail from previous suggestions that the 
gas consumption timescale was important for elliptical FP 
properties (e.g. Bekkil ll998h . Other theoretical studies of el- 
liptical galaxy formation have suggested that gas dissipation 
was necessary for sett ing or ma intaining elliptical galaxy scal- 
ing laws, notably Ciotti & van Albada (2001) who conclude 
from dissipationless mergers that significant dissipation is re- 
quired to simultaneously satisfy both the Mbh-o" and FP rela- 
tion. Combined with the results of D i Matteo et al..(k200 5) and 
Robertson et al. (2005a), this work demonstrates that highly 
dissipational mergers can generate both the FP and Mbh-c 
relations with tight scatter. Our analysis has also provided a 
detailed explanation of how dissipation systematically varies 
with galaxy mass and to what extent it induces FP tilt. 

Furthermore, our dissipationless simulations are consistent 
with the emerging picture that at least one generation of 
spheroidal mergers do not strongly invalidate the scaling laws 
of their elhptical galaxy prog enitors (Capelato et al. 1995t 
iDant as et al1l2003MNipoti et aL .2003: Bovlan-Kolchin et all 
2005). The recent observations suggesting that elliptical 
galaxies undergo 0.5 - 1 such mergers at redshifts z < 1 
(Belletal. 2005; van Dokkum 2005) may require such be- 
havior as the tight local FP relation observationally limits the 
possible importance of scenarios that substantially increase 
the FP scatter We also note here that the results from our 
simulations of the Mbh-c relation produced by re-mergers 
of ellipticals are in rough agreement with previous studies 
(e.g. Nipoti et al. 2003). Moreover, the fact that re-merger 
remnants do not increase in velocity dispersion dramatically 
while their black hole masses double suggests that if black 
holes are overly massive relative to their host bulges at high 
redshifts (e.g. Walter et al. 2004; Peng et al. 2005), then major 
dissipationless mergers of spheroids will only exacerba te their 
discrepancy from the loca l MRH -g relation ( Gebhar dt et alJ 
llOOO; Ferrarese & Merritti l200(i |Tremaine et al. 200|). If 
the structural properties of quiescent galaxies instead in- 
crease the central velocity dispersions at a given black hole 
mass at high redshifts, major dissipationless mergers of 
spheroids will decrease the scatter in the M^n-a relation lo- 
cally (Robertson et al. 2005b). Lastly, if the frequency of ma- 
jor dissipationless mergers increases with galaxy mass then 
the high-ma ss end of t he Mbh-c relation may steepen relative 
to the Tremaine et al. (2002) best-fit rel ation and a ppear as a 
mass-dependent scatter in the data (Robertson et al. 2005b). 

Independent evidence from simulations that dissipation is 
important for the properties of elliptical galaxies is mount- 
ing. Using dissipational simulations of disk galaxy mergers, 
including a subset of the simulations analyzed in this work, 
£qx et al. (2005 b) demonstrate that progenitor gas fractions 
of order /gas ~ 0.3 or higher are necessary to reproduce the 
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kinematic structure of elliptical galaxies as a function of their 
stellar mass. Notably, the vja ratio of rotation to dispersion 
of remnants from /gas = 0.4 mergers better reproduces the ob- 
served distribution of elliptical kinematic s than do purely dis - 
sipationless disk galaxy mergers (e.g. Naab & Burkert 2003J). 
Our analysis demonstrates that the Cox et al. (2005b ) rem- 
nants also lie on a tight FP relation. 

While recent observational determinations of the fun- 
damental plane in galaxies at redshifts z ^ 1 suggest 
that the MjL offset of the FP changes with redshift 
owing to passive evolution of the stellar populations, 
there is still disagreement reg arding the r e dshift- 

evolution of the FP tilt (e. g. iGe bhardt et al.' '200^ 

Ivan Dokkum & Stan ford! 120031 Ivand er Wei et al. 20ol 
l^eu et al. 2005; di Serego Alighieri et al. 2005). The results 
of our simulations of disk galaxy mergers appropriate for 
higher redshifts suggest that if the MjL tilt induced by 
stellar population effects remain fixed as a function of galaxy 
mass with redshift (i.e. evolution in the stellar populations 
only affects the FP zero-point), then the observed FP tilt 
will remained fixed out to high redshifts. The implication 
from stellar population modeling that elliptical galaxies 
likely formed early (z > 1-2) then may not provide a 
tight constraint on the merger hypothesis as a theory for 
elliptical galaxy formation since gas-rich disk galaxy models 
appropriate for a large range in redshift produce the same 
stellar-mass FP relation. As methods for including a cosmo- 
logical accounting of elliptical galaxy formation times and 
stellar population color evolution into structu re formation 
theory improve (see, e.g. 'Ho pkins et al ]'2005b', for a related 
methodology), theoretical models will have more to say in 
the future on this topic. 

As a final note, some recent observations have been in- 
terpreted as evidence that the FP tilt is consistent with a 
purely stellar-population induced MjL variation with galaxy 
mass (e.g. fCappellari et al. 2005). The success of Cox et aj 
(|2005b) at reproducing the kinematic properties of elliptical 
galaxies with gas-rich disk galaxy merger simulations (includ- 
ing some simulations analyzed in this paper) argues strongly 
that dissipation is important for understanding kinematical 
mass trends. Our simulations suggest that a corollary of the 
dissipational origin of kinematic structure in ellipticals is FP 
tilt, which at least partially owes to central Mtotai/A^* varia- 
tions with galaxy mass. We note that this tilt may be in addi- 
tion to both kinematic nonhomology-related and photometric 
nonhomology-related tilt, and a full characterization of these 
effects is planned for future work. 

7. SUMMARY 

We determine the Fundamental Plane and R^-M^, relations 
produced by the merging of galaxies appropriate for redshifts 
z = 0-6. We demonstrate that gas dissipation induces tilt into 
the Fundamental Plane, causing a deviation from the theoret- 
ical scalings derived from the virial theorem for homologous 
systems. Furthermore, the recovery of the R^-M^, relation for 
early-type galaxies similarly requires gas-rich progenitors to 
steepen the relation relative to the R^-M^, relation produced 
by the dissipationless merging of disk galaxies. These simu- 
lations argue that the merger hypothesis for the formation of 
locally-observed elliptical galaxies is successful if disk galax- 
ies merge when dissipational effects are sufficiently impor- 
tant, requiring gas fractions /gas > 0.3 in disk galaxies at red- 
shifts where ellipticals are primarily formed (typically z > 1). 
We provide a detailed summary of these results below. 



• The Fundamental Plane (FP) relation produced by the 
merging of dissipationless disk galaxy models lies 
nearly parallel to the plane defined by the virial relation 
for homologous systems. Increasing the angular mo- 
mentum of the orbit by lengthening the pericentric pas- 
sage distance of the orbit produces an offset in the FP by 
increasing the effective radius of the remnants. Scaling 
the disk galaxy models for redshifts z = - 6 produces 
the same FP relation. The effective radius R^ - stel- 
lar mass relation produced in these mergers is shal- 
lower than that measured for massive galaxies in the 
Sloan Digital Sky Survey (Shen et al. 2003). Wide or- 
bits and progenitors appropriate for low-redshifts pro- 
duce larger remnants than do high- redshift progeni- 
tors or progenitors with bulges, but all dissipationless 
mergers produce shallow R^-M^, relations. We also use 
higher resolution simulations to demonstrate that these 
results are insensitive to our numerical resolution. 

• The merging of gas-rich disk galaxy models with cool- 
ing, star formation, and supernova feedback produces 
a Fundamental Plane (FP) scaHng {R^ oc cr'-5'V;° '*°) 
similar to the observed infrared FP (R^ oc cr'^^/"'^ ™; 
Ifahre et al. 1998b) and is almost independent of the 
redshift scalings of the progenitor systems. The ef- 
fective radius R^. - stellar mass relation produced 
by these mergers is roughly parallel to that measured 
for massive galax ies in the Sloan Digital Sky Survey 
dShen et a l1 l20?)3l) . but the chosen prograde-prograde 
coplanar orbit generates remnants that lie somewhat be- 
low this relation. Higher-redshift progenitors produce 
smaller remnants. 

• The Fundamental Plane (FP) relation produced by the 
merging of gas-rich disk galaxies with cooling, star 
formation, supernova feedback, and a prescription for 
feedback from accreting supermassive black holes ex- 
hibits a scaling (/?e (t''''^/J° **^) similar to the observed 
infrared FP (R^ oc cr' "4"" ™; Pahr e et all 199 8b) and is 
nearly coincident with the FP produced by similar sim- 
ulations without black holes. The FP relation is roughly 
independent of the redshift scalings of the progenitor 
systems and the location of remnants within the FP is 
fairly insensitive to a large variety of disk orientations 
and orbital configurations as changes in the effective 
radius are compensated by changes in the velocity dis- 
persion and surface mass density. The effective radius 
Rf. - stellar mass relation produced by the same sim- 
ulations produces a Rs-M^, relation roughly parallel to 
that measured for massive galaxies in the Sloan Digi- 
tal Sky Survey (Shen et al. 2003), with an offset toward 
smaller remnants. Varying the system angular momen- 
tum through combinations of the initial disk orientation 
and pericentric passage distance for a single pair of pro- 
genitor models produces a spread in the remnant effec- 
tive radius, and a proper accounting of cosmologic alor- 
bits would decrease the discrepancy with the She n et alJ 
(2003) normalization. The Ri^-Mi, relation shifts to- 
ward smaller remnants as the redshift of the progenitor 
systems increases. 

• We estimate that « 40-100% of the Fundamental Plane 
(FP) tilt produced in the simulations owes to an increase 
in stellar mass at high phase-space densities and bind- 



20 



Robertson et al. 



ing energies as gas converts into stars in the central- 
most regions of the remnants. This increase in phase- 
space density corresponds to a decrease in the total- 
to-stellar mass ratio Mtotai/A/* in the central regions 
of simulated elliptical galaxies. Furthermore, the de- 
crease in Mtotai/A/* is a function of galaxy mass, with a 
stronger decrease in smaller systems, thereby inducing 
a FP tilt. We find that the gas-rich mergers (/gas > 0.3) 
required to steepen the R^ -M^, relation to become par- 
allel with the observed Shen et aP J2003i) relation also 
produce significant FP tilt. 

• The mass-dependent impact of dissipation on the cen- 
tral Mtotai/A/* of galaxies is tied to scale-dependent pro- 
cesses in the merging systems. In massive halos, a 
substantial fraction of the available gas is collisionally 
heated during the first pericentric passage and driven 
into the halo. This heated gas cannot cool efficiently, 
thereby causing massive systems to act more like dis- 
sipationless systems than low mass merging galaxies. 
Moreover, the dynamic variation of the consumption 
timescale in low mass galaxies allows for the rapid con- 
version of gas into stars during only the most violent in- 
teractions of the merger while conserving the gas con- 
tent of the system in more quiescent phases and allow- 
ing for gas dissipation to have an increased impact on 
the central properties of the remnants. 

• The Fundamental Plane (FP) relation produced by the 
merging of spheroidal Hernauist ( 1990) galaxy mod- 
els with dark matter halos remains similar to the FP 
delineated by their progenitors, in agreement with pre- 
vious results (e.g. Bovlan-Kolchin et al. 2005). Addi- 
tionally, the re-merging of the remnants from gas-rich 
disk galaxy mergers also roughly maintains the FP of 
their progenitors. The effective radius - stellar mass 
Mi, relations produced by the merging of spheroidal 
galaxy models produces a slightly shallower R^ - 
relation than that delineated by the progenitor models, 
and although the remnants nearly double in effective 
radius through equal-mass mergers, a single generation 
of remnant re-merging will likely induce scatter but not 
destroy the Rf.-Mi, correlation. 

• The re-merging of ellipticals formed from the merg- 
ing of disk galaxies will likely induce additional scat- 
ter in the local black hole mass (Mbh) - stellar ve- 
locity dispersion ia) relation, but should not destroy 
the correlation after a single generation of remnant re- 
merging. In agree ment with previous simulations (e.g. 
tNipoti et al.l2003l) . we find that the velocity dispersions 
of the re-merger remnants do not increase substan- 
tially relative to their progenitors during binary merg- 



ers. However, the typical number of 0.5 - 1 dissipa- 
tionless mergers below z = 1 recently inferred frorn 
observations (e.g. ' Bell et al.ll2005t Ivan Dokkumll2005t) 
should not dramatically change the Mbh-c scaUngs im- 
printed at higher redshifts. If the velocity dispersion 
in remnants for a given black hole mass is larger at 
high redshifts than that observed locally, as inferred 
from simulations of the generation of the Mbh-c re- 
lation fro m the merger of disk galaxies at redshifts 
z = 0-6 (Robertson et al. 2005b), the dissipationless 
re-merging of remnants would decrease the scatter in 
the local Mbh-c relationship under the condition that 
the strength of black hole feedback did not vary sig- 
nificantly with redshift. If black holes are substantially 
more massive than their host bulges at high redshifts 
(e.g. Walter et al. 2004; Peng et al. 2005), major dissi- 
pationless mergers of spheroids at low redshifts would 
exacerbate the discrepancy of those galaxies with the 
local Mbh-c relation. An increasing frequency of 
major dissipationless mergers with galaxy mass may 
steepen the high-mass Mm-a relation relat ive to the 
local best-fit trend (e.g. lTremaine et al.l2002h . 



The success of the merger hypothesis for elliptical galaxy 
formation seems intimately connected with the physical pro- 
cess of dissipation. Our results demonstrate that if disk galaxy 
mergers are strongly dissipational, tilt is generated in Funda- 
mental Plane and the R^-M^, relation steepens to approach the 
observed scaling. Combined with compelling evidence from 
the kinematic s tructure of simulated elliptical galaxies (e.g. 
ICox et al. 2005b). we conclude from our results that the for- 
mation of elliptical galaxies most likely occurred in gas-rich 
mergers of disk galaxies. Reconciling observational determi- 
nations of the Fundamental Plane with the merger hypothesis 
will require the abandonment of simplistic assumptions for 
the relative stellar and dark matter mass distributions at the 
centers of elliptical galaxies for the interpretation of observa- 
tions and require more sophisticated modeling of cosmologi- 
cal structure formation, including a full cosmologically moti- 
vated distribution of elliptical stellar colors and metallicities. 
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